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Abstract 



We have studied the star-formation and AGN activity of massive galaxies in the redshift range z = 0.4 — 2, 
which are detected in a deep survey field using the AKARI InfraRed (IR) astronomical satellite and 
Subaru telescope toward the North Ecliptic Pole (NEP). The AKARI/IRC Mid- InfraRed (MIR) muhiband 
photometry is used to trace their star-forming activities with the Polycyclic- Aromatic Hydrocarbon (PAH) 
■ emissions, which is also used to distinguish star-forming populations from AGN dominated ones and to 

' estimate the Star Formation Rate (SFR) derived from their total emitting IR (TIR) luminosities. In 

^ . combination with analyses of their stellar components, we have studied the MIR SEE) features of star- 

forming and AGN-harboring galaxies, which we summarize below: 1) The rest-frame 7.7-/im and 5-/im 
luminosities are good tracers of star-forming and AGN activities from their PAH and dusty tori emissions, 
respectively. 2) For dusty star-forming galaxies without AGN, their SFR shows a correlation that is 
nearly proportional to their stellar mass and their specific SFR (sSFR) per unit stellar mass increases with 
redshift. Extinctions estimated from their TIR luminosities are larger than those from their optical SED 
fittings, which may be caused by geometric variations of dust in them. 3) Even for dusty star-forming 
galaxies with AGN, SFRs can be derived from their TIR luminosities with subtraction of the obscured 
AGN contribution, which indicates that their SFRs were possibly quenched around z ~ 0.8 compared with 
those without AGN. 4) The AGN activity from their rest-frame 5-/xm luminosity suggests that their Super 
Massive Black Holes (SMBH) could already have grown to ~ 3 x IO^Mq in most massive galaxies with 
IO^^Mq at z > 1.2, and the mass relation between SMBHs and their host galaxies has already become 
established by z ~ 1 — 2. 

Key words: galaxies: evolution — galaxies: high-redshift — infrared: galaxies 
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1. Introduction 

A major issue in observational cosmology is to recon- 
struct the cosmic evolution of star formation and the de- 
tailed growth history of the present-day galaxies. Even 
though it is not easy to reconstruct the whole path of 
the stellar mass for galaxies as a function of the present- 
day galaxy stellar mass M^.q and the cosmic time t as 
M*(M*;o,t), it is essential to obtain knowledge of their 
Star Formation Rate (SFR) as a function SFR{M^,t). 

Recent near-infrared (NIR) observations suggest that 
most of the stellar mass is contained in massive galaxies, 
and that this situation had already become established at 
z ~ 1 . This appears to have been followed by a rapid evolu- 
tion of the stellar mass density between z = 1 — 3 (Papovich 
et al. 2001; Dickinson et al. 2003; Fontana et al. 2003; 
Franx et al. 2003; Drory et al. 2004; Fontana ct al. 2004; 
Glazebrook et al. 2004; van Dokkum ct al. 2004; Labbc 
et al. 2005; Fontana ct al. 2006; Papovich 2006; Yan ct al. 
2006; Webb et al. 2006; Kriek et al. 2006; Pozzetti et al. 
2007; Arnouts et al. 2007; Vcrgani et al. 2008). These re- 
cent studies imply that most of the massive galaxies must 
already have formed at z ~ 1, supporting the so-called 
downsizing scenario (Cowie et al. 1997). In this scenario, 
star formation shifts from higher to lower mass galaxies as 
a function of cosmic time t, although the stellar mass in 
massive galaxies may not directly follow the growth his- 
tory of their Dark Matter Halos (DMHs) which build up 
hierarchically from less massive systems to more massive 
ones at a later time. However, several studies of massive 
galaxies at z < 1 (Brown et al. 2007; Faber et al. 2007; 
Brown et al. 2008) report the evolution from predomi- 
nantly blue to red galaxies, suggesting that the number 
and stellar mass of galaxies in the red-sequence may con- 
tinue to grow, even after z = 1 (Bell et al. 2004; Zucca 
et al. 2006), in agreement with the hierarchical merging 
scenario. Recently, a mass-dependent quenching of star 
formation has been proposed to explain these rapid in- 
crease of the number density of the massive galaxies rela- 
tive to lower-mass galaxies at z = 1 — 2 (Peng et al. 2010; 
Kajisawa 2011; Brammer et al. 2011). 

The evolutionary path of the M*(M*;o,i) up to the 
present mass Af*;o can be also reconstructed from knowl- 
edge of the SFR{AI^,, t). Recent studies have con- 
firmed the features for rapid SFR evolution up to z ~ 2 
(Brinchmann & Ellis 2000; Fontana et al. 2003; Perez- 
Gonzalez et al. 2005; Fculncr et al. 2005a, b; Papovich 
2006; Iglesias-Pramo et al. 2007; Zheng ct al. 2007; Cowie 
& Barger 2008; Damen et al. 2009; Dunne et al. 2009; 
Santini et al. 2009; Pannella et al. 2009). However, we 
have not found a general consensus concerning the stel- 
lar mass dependence for specific SFRs: sSFR = SFR/M^ 
(the SFR per unit stellar mass). Most studies based on 
UV derived SFR's have found evidence to support an anti- 
correlation with stellar mass, not only local galaxies, but 
also the distant Universe, supporting the downsizing sce- 
nario. It is not easy, however, to accept that UV estimates 
are the only way to estimate the SFR that will capture all 
the emission originally from new born stars in a galaxy. 



which can be independently estimated from rest-frame lu- 
minosities in the UV, IR, Submillimeter (Submm), and 
radio wavelength regions by introducing appropriate cali- 
bration factors. Although the UV luminosity is used as a 
proxy of the SFR, its reliability often depends on the use of 
large, and rather uncertain corrections for dust extinction. 
Most of the short wavelength photons from star-forming 
galaxies are reprocessed through absorption and scatter- 
ing with dust in their star-forming regions, and then re- 
emittcd in the IR- Submm wavelength regions. Thus, the 
mid/far-IR and their submillimeter emissions become an 
important tracer of the total star forming activity in high- 
z galaxies. 

The total IR (TIR) luminosity Ljb. integrated across 
the long wavelength spectral region is a good proxy 
for the SFR, which is relatively independent of the 
dust obscuration (Kennicutt 1998). Luminous InfraRed 
Galaxies (LIRGs) and Ultra Luminous InfraRed Galaxies 
(ULIRGs) up to z = 1 - 2, detected at the Mid-InfraRcd 
(MIR) wavelength with the IR observing satellites as the 
ISO (Mann et al. 2002; Elbaz et al. 2002; Pozzi et al. 2004; 
Dennefeld et al. 2005; Taylor et al. 2005; MarciUac et al. 
2006b,a), the AKARI (Takagi et al. 2007, 2010), and the 
Spitzer (Egami ct al. 2004; Shapley et al. 2005; Caputi 
et al. 2006; Choi ct al. 2006; Marcillac et al. 2006a,b; 
Papovich 2006; Papovich et al. 2006; Rcddy et al. 2006; 
Elbaz et al. 2007; Papovich et al. 2007; Rocca-Volmcrangc 
et al. 2007; Shim et al. 2007; Teplitz et al. 2007; Brand 
et al. 2008; MagneUi et al. 2008; Dasyra et al. 2009; 
Hernan-Caballero et al. 2009; Murphy et al. 2009; Santini 
et al. 2009; Strazzullo et al. 2010), have shown that the 
SFR rapidly increases with the redshift, as they contribute 
~ 70% of the cosmic infrared luminosity density (Floc'h 
et al. 2005). The MIR emission from LIRGs/ULIRGs 
is evidence of dusty starburst or AGN surrounded by a 
dusty torus. Recently, the LIRGs/ULIRGs are expected 
to suppress star formation at z = 1 — 2 with feedback ac- 
tivity from starbursts or AGNs. Understanding of the 
LIRGs/ULIRGs at z = 1 — 2 can be a key to solve the 
cosmic history of star formation and AGN activities. The 
SFR for galaxies detected with the Spitzer is frequently 
estimated with the TIR luminosity estimated from the 24 
/im flux (Daddi et al. 2004; Santini et al. 2009; Pannella 
et al. 2009), which is nearly proportional to their stel- 
lar mass up to z ~ 2. These recent results with the 
Spitzer are also basically consistent with recent results 
from stacking analysis for their radio emissions (Dunne 
et al. 2009; Pannella et al. 2009) as another SFR estima- 
tion immune to dust extinctions. These two radio stack- 
ing analyses, however, show a discrepancy in their details 
as Dunne et al. (2009) report sSFR with a negative cor- 
relation to stellar mass while Pannella et al. (2009) re- 
port sSFR with stellar mass independence. Thus, there 
may be various unknown factors in the SFR estimations. 
The TIR luminosity Ljr of LIRGs/ULIRGs is generally 
derived from the fitting of the observed Spectral Energy 
Distribution (SED) with synthetic templates. In the SED 
fitting, we could sometimes use only a few bands in the 
IR wavelength without longer 100 — 1000 /xm wavelengths, 
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which may cause large uncertainties in the Lm estima- 
tion. Obscured AGNs are another notable problem in the 
Lju estimation. It means that a sample selection without 
significant AGN contamination is also important to ob- 
tain an accurate SFRjji. Thus, systematically designed 
surveys are required to construct samples for distant star- 
forming galaxies. 

Optical surveys have been successful in selecting numer- 
ous Lyman Break Galaxies (LBGs) as star forming galax- 
ies beyond z ^ 2 (Steidel et al. 1996a, b, 1999; Lehnert 
& Bremer 2003; Iwata et al. 2003; Ouchi et al. 2003, 
2004; Dickinson ct al. 2004; Sawicki & Thompson 2006a,b) 
and even up to z ~ 6 (Yan et al. 2004; Stanway et al. 
2003; Bunker et al. 2004; Bouwens et al. 2004; Shimasaku 
et al. 2005). They may, however, fail to detect substan- 
tial numbers of massive galaxies at 2 < 3, in which the 
Lyman break photometric technique does not work well 
with ground-based telescopes because of telluric absorp- 
tion. Furthermore, even if galaxies are undergoing inten- 
sive star formation, the internal dust extinction of their 
UV radiation may reduce the prominence of their Lyman 
break features, to appear as LIRGs/ULIRGs in the local 
Universe. However, alternative photometric techniques 
can be used to select massive high-z galaxies through their 
Balmer/4000 A break or the 1.6 /im bump features in their 
Spectral Energy Distributions (SEDs), as evidence of the 
presence of an evolved stellar population. We will desig- 
nate the former and the latter as 'Balmer' Break Galaxies 
(BBGs) and as IR 'Bump' Galaxies (IRBGs), respectively. 
It is known that both star-forming and passive BBGs in 
the redshift desert around z = 2 can be selected with a 
simple and sophisticated two-color selection based on B,z, 
and K band photometry (Daddi et al. 2004). This three 
filter optical-NIR technique was designed to detect the 
Balmer/4000 A break, nearly free from the dust extinction 
which caused some difficulties for previous 'red color' tech- 
niques relying on a color diagnostic of {R — K)vega > 5 — 6, 
corresponding to {R — K)ab > 3.35 — 4.35 (Elston et al. 
1988; Thompson et al. 1999; Daddi et al. 2000; Firth et al. 
2002; Roche et al. 2002, 2003) or for (J ~ K)vega > 2.3 
corresponding to (J — K)ab > 1-3 ^ (Pozzetti et al. 2003; 
Franx et al. 2003). With the ncar-IR (NIR) observations 
from the Spitzer and from the AKARI, it has become pos- 
sible to select IRBGs at z > 1 with NIR colors at > 3 fim 
of the 1.6 ^m bump (Simpson & Eisenhardt 1999; Sawicki 
2002; Berta et al. 2007). These optical-NIR observations, 
detecting the Balmcr/4000 A break and the 1.6 /^m bump, 
have been also useful to trace how stellar mass is corre- 
lated with the absolute NIR luminosity. 

Although optical-NIR multi-band photometry has the 
potential to select a sample for star-forming massive 
galaxies, it may not be completely free from the AGN 
contamination. MIR multi-band photometry tracing PAH 
emissions at 6.2, 7.7 and 11.3 /im and Silicate (Si) ab- 
sorption at 10 fim (rest-frame) from dusty star-forming 
regions, can be used to discriminate between the star- 
forming activity in LIRGs/ULIRGs and that of AGN 

^ Jab = Jvega + 0-89, KSAB = Ksvega + 1-84 



dominated galaxies. Thus, systematic multi- wavelength 
observations that include observations at > 2 /tm are 
important for studies the star-forming activity of mas- 
sive high-z galaxies, and in looking back at the evolu- 
tion of their SFRs. In order to connect the history of 
stellar mass growth in massive galaxies traced in optical- 
NIR wavelength with ground based telescopes and that of 
star formation traced in MIR wavelength observed out- 
side of the atmosphere, we undertook multi-wavelength 
surveys of deep and wide areas with ~ 0.38 and ~ 5.8 
deg^, respectively, in a field close to the North Ecliptic 
Pole (NEP), which were extensively observed at the IR 
wavelength with the first Japanese infrared astronomical 
sateUite AKARI. 

In this paper, we attempt to reconstruct the cosmic his- 
tory of dusty massive galaxies up to z ~ 2, selected esti- 
mated from optical-IR images in the NEP survey field. We 
report on: 1) photometric redshift Zphot, stellar mass M*, 
and SFRuv;cor, which arc estimated with their optical- 
NIR SED fittings, 2) the distinction between dusty star- 
burst and AGN dominated galaxies using the MIR SED 
feature, 3) TIR luminosities of starburst dominated galax- 
ies with the subtraction of the obscured AGN contribu- 
tion, from which we can derive S F Rj fi^uv complemen- 
tary to the SFRuv.cor; and 4) IR luminosities of the dusty 
tori of AGN dominated galaxies, which show evolutionary 
trends in the connection between AGN and star forming 
activities. 

In section 2, we summarize the data we used; in sec- 
tion 3, we present the photometry for z'-band detected 
galaxies in our multi-wavelength dataset; in section 5, we 
also report AKARI/IRC detected IRBGs classified using 
the rcdshiftcd 1.6-/im IR bump in the NIR bands, which 
are divided into three rcdshiftcd subgroups; in section 6, 
we present AKARI/IRC MIR detected galaxies which cor- 
respond to LIRGs/ULIRGs and classify them into dusty 
starbursts and AGN from their MIR SEDs with/without 
PAH features in the MIR bands; in section 7, we derive 
the TIR luminosities for the MIR detected galaxies; in 
section 8, we estimate their star formation rates and ex- 
tinctions; in section 9, we reconstruct their intrinsic stel- 
lar populations with extinction correction; in section 10 
we also study the AGN activities in the MIR detected 
galaxies; in section 11, we discuss some expectations for 
on-going and future observations; finally in section 12, we 
summarize the conclusions of our study. 

In this paper we adopt Hq = 70 km s~^ Mpc~^, = 
0.3, = 0.7 throughout. All magnitudes are represented 
in the AB system (e.g. Oke 1974). 

2. The data set 

We have used mainly the data sets obtained in the 
AKARI Deep imaging survey with all of the available 
bands of the Infra Red Camera (IRC) on the AKARI 
(Onaka et al. 2007) near the North Ecfiptic Pole (NEP). 
Wc define the N2, Ni, iV4, ST, SQW, 5*11, L15, LISW , 
and L24 bands as having central wavelengths at 2, 3, 4, 7, 
9, 11, 15, 18, and 24 /xm respectively. These are reported 
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Table 2. Summary of limiting magnitudes and photometric 
corrections. 
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in this paper as the IR photometric AB magnitude of an 
object observed in the various IRC bands. Details of the 
data reduction of the NEP Deep and Wide surveys have 
been presented by Wada et al. (2UU8) and Lee et al. (2007) 
(see also Hwang et al. (2007) and Ko et al. 2012), respec- 
tively. We have used the same IRC scientific data for the 
AKARI NEP Deep field as described in detail by Wada 
et al. (2008), which include mosaicing them together and 
making pointing and distortion corrections, as well as cor- 
rection to mitigate detector artifacts. 

The AKARI NEP Deep survey area has been cov- 
ered by observations in the near-UV u* band with the 
CFHT/MegaCam(M-Cam) (Boulade et al. 2003); the op- 
tical B,V,R^i', and z' bands with the Subaru/Suprimc- 
Cam (S-Cam) (Miyazaki et al. 2002); and in the NIR J 
and Ks bands with the KPNO(2.1m)/Flamingos(FLMG) 
(Elston et al. 20U6). Details of the Subaru/S-Cam, 
CFHT/M-Cam, and KPNO/FLMG data reduction can 
be found in Imai et al. (2007), Takagi et al. (2012), 
and Ishigaki et al. (2012). In this study, we have se- 
lected ~ 56000 galaxies with z' < 26.4 and S/N> 3. For 
these sources, we estimate their multi-wavelength photo- 
metric characteristics, adding ground-based near-UV, op- 
tical, and NIR imaging data to the AKARI/IRC data as 
described in the following Section. For simplification, we 
use u,i,z, and K as u*,i',z' and K^, respectively for rep- 
resenting colors derived from the respective filters. 

3. Photometry 

3.1. Optical and NIR photometry for ground-based data 

Object detection and photometry were made using 
SExtractor 2.3.2 (Bertin & Arnouts 1996). The z' S- 
Cam image was used to detect optical objects, defined as a 
source with a 5 pixel connection above the 2 a noise level. 
We derived the Full Width Half Maximum (FWHM) of 
the Point Spread Function (PSF) from the photometric 
images of stellar objects as point sources (see their se- 
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" m,iim. is the hmiting magnitude. 

Arria.c. is the aperture correction in magnitude. 
h.mg,s, is the magnitude correction for galactic 
the field. 

Arris. o. is the magnitude correction for systemat: 
the SED calibration for the stellar objects. 



extmction m 
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lection scheme in subsection 3.4), which are summarized 
in table 1. In the process to make a catalogue, the rea- 
son why the z'-band was chosen as the detection band 
is the following: 1) the FWHM of the PSF 1" in the z' 
image with the S-Cam is good enough to prevent the con- 
fusion, 2) the z' photometry is one of the essential bands 
in the star-galaxy selection with the BzK color as shown 
in subsection 3.4, 3) most of the galaxies of interest as de- 
tected by the AKARI; N3 Red Galaxies (N3Rs) and MIR 
bright N3Rs (MbN3Rs), and optical- NIR selected galax- 
ies as BBGs are clearly identified in the z' image. See 
the details about the N3Rs, the MbN3Rs, and BBGs in 
subsections 5.1 and 6.1, and in appendix 1), respectively. 
The catalogue of detected objects to z' < 26.4 with a 2" 
diameter aperture was generated over the whole 23' x 32' 
area of the S-Cam z' image. 

After selecting ^ 56000 z' sources in the area observed 
with the FLMG, from the ^ 66000 sources detected in 
the S-Cam z' images, we performed aperture photometry 
in the near-UV u*, optical B,V,R,i\z', and NIR J, Kg 
images, fixed at the positions of the z'-detected sources. 
The photometry on the B,V, R,i' , and z' images was 
performed with a small 2" diameter aperture with the 
SExtractor. The photometry on the u* and the J and 
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Ks images was performed with larger 3" and 4" diame- 
ter apertures with the photo of the IRAF task since these 
image quahties in off-center regions are not good due to 
the aberrations compared with those of the S-Cam. The 
larger aperture photometry reduces the uncertainties in 
the photometry and the alignment errors in the World 
Coordinate System (WCS) from these aberrations. The 
errors of the ground-based optical and NIR photometry 
come mainly from the sky background and its fluctuation. 
We have estimated them from a Monte Carlo simulation 
with the same method as that of Ishigaki et al. (2012), 
in which we measured the counts in randomly distributed 
apertures and fit the negative part of the count histogram 
with a Gaussian profile. Thus, the Scr limiting magni- 
tudes TOiim were derived from the background estimation, 
as summarized in table 2. 

The final integrated magnitudes for all filters were ob- 
tained after applying an aperture correction AttIq c. based 
on empirical PSFs, which are constructed from the pho- 
tometry for stars in each image band. All of the optical 
and NIR magnitudes were further corrected for a galactic 
extinction of E{B — V)= 0.041, based on observations of 
the survey area taken from Schlegel et al. (1998) and us- 
ing the empirical selective extinction function of Cardelli 
et al. (1989). 

3.2. Infrared photometry for AKARI/IRC images 

As shown above, we have produced a multiband merged 
catalogue from the ground-based telescope images for the 
z'-detected objects in the field. As shown in table 2, the 
depth of A^2 images obtained with the AKARI/IRC is 
comparable to that of Kg images with the KPNO/FLMG, 
although the images of all other bands at longer wave- 
lengths are shallower than those of z'. In fact, from 5819 
sources extracted with the SExtractor in the N3 band 
with a 5 pixel connection above the 2 a noise level, 92% of 
them can be identified as 5332 counterparts of z'-detected 
sources without confusion. Thus, we have performed the 
photometry for all of the AKARI/IRC bands on the z'- 
detected sources. 

The IRC photometry on the z'-detectcd sources was 
performed with 8" and 10" diameter apertures for NIR 
{N2,N3, and A^4) and MIR {S7,S9W,S11,L15,L18W, 
and L2A) bands, respectively, with the IRAF task photo. 

The final integrated magnitudes in the IRC photometry 
were obtained after applying an aperture correction based 
on the growth curves with empirical PSFs, constructed 
from photometry of stars in each band image, which is 
basically similar to the scheme descirbed in Takagi et al. 
(2012). The aperture corrections for the IRC are sum- 
marized in table 2, along with those of the ground-based 
near-UV, optical, and NIR magnitudes. 

In order to reduce the source confusion due to the fact 
that the IRC PSFs are notably worse than those in the 
ground-based observed images, we have performed IRC 
photometry with several different apertures in each IRC 
band, and compared these magnitudes, excluding objects 
with a large magnitude difference from the IRC counter- 
parts in the z'-detected catalogue. In the photometry in 



N2,N3, and A^4 bands, we have optimized the constraint 
for the IRC/NIR counterparts as Af3(6") -iV3(8") < 0.51 
between 6" and 8" diameter photometry in the N3 band, 
adding another constraint of Ks — A^2 < 1.0 at 2 iim. 

For MIR sources detected with SN3 at least in one 
MIR band, there were many more sources with confu- 
sion compared with the NIR sources. We could select ~ 
800 MIR bright N3 Red galaxies (MbN3Rs) and - 900 
MIR marginally-detected N3 Red galaxies (MmN3Rs) as 
shown in section 6. However, we have taken confirmed 
samples of ~ 600 MbN3Rs and 600 MmN3Rs after ex- 
cluding NIR counterparts of the MIR sources having pairs 
within 6", in order to correctly select the candidates of 
LIRGs at z > 0.4 as discussed in section 6. Roughly two 
thirds of these MbN3Rs and MmN3Rs were also iden- 
tified as sources listed in another MIR source catalogue 
reported by Takagi et al. (2012), which are almost the 
same after taking the difference in their detection lim- 
its into consideration, as the former was selected with 
SN > 3 at least in one MIR band while the latter was 
with SN > 5. Concerning these NIR and MIR detected 
galaxies, we could confirm with our eyes that it can work 
to mostly exclude confusion from neighborhood sources 
within a circleregion corresponding to the IRC PSF ra- 
dius around each z' source. 

3.3. Photometric check and systematic offsets 

In order to further check the photometric zero points, 
we have performed SED fitting with u* BVRi' z' JKg and 
AKARI N2,N3, and iV4 photometric data for the steUar 
objects. It is well known that galactic stars can be clearly 
selected using a two-color criterion (Daddi ct al. 2004; 
Kong et al. 2006; Lane et al. 2007; Quadri et al. 2007; 
Blanc et al. 2008): 

BzK, = {z- K)ab - 0.3(B - z)ab < -0.3 , (1) 

where we have taken —0.3 as a parameter in the right 
side, which defines a line separating stars and galaxies 
that is slightly larger than Daddi's original value of —0.5, 
to reduce the contamination of stars in our galaxy sam- 
ple. We have used the BzK color criterion to select the 
photometric calibrators as stellar objects with additional 
criteria for the magnitudes z' > 19 and stellarity parameter 
CLASS_STAR > 0.95 in SExtractor (see also subsection 
3.4). We have also excluded IRC/NIR confused image 
sources by visual examination. 

We have adopted SED templates from the stellar library 
of Bruzual-Persson-Gunn-Strykcr (BPGS) Atlas^, based 
on Gunn & Stryker (1983), which is similar to the Pickles 
(1998) stellar spectra'^, and Lejeune et al. (1997). In order 
to include the 7V4 magnitudes in our fitting of SEDs to 
the stellar objects with the BPGS Atlas, we have extended 
the library SEDs beyond wavelengths of A > 2.5 /im with a 
Rayleigh- Jeans law. The best fit is sought by means of 
minimization, among all templates in the stellar library 
after varying their normalization. We have estimated the 

http:/ /www. stsci.edu/hst/obsorvatory/cdbs/bpgs. html 
http:/ /www. ifa.hawaii.edu/uscrs/pickles/AJP/hiHb. html 
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offsets in each of the bands, which are seen as systematic 
discrepancies between the original observed magnitudes 
and the best-fitted library magnitudes. These systematic 
offsets are also summarized in table 2. Since the derived 
offsets are less than 0.2, basically we can confirm the ac- 
curacy of the photometric zero points in all the bands. In 
the comparison between the spectroscopic and photomet- 
ric redshifts, the result for the systematic offset is slightly 
better than that without this. Thus, we use the offset. 

3.4- Star-galaxy separation 

As discussed in subsection 3.3, the colors estimated with 
the BPGS star atlas overlap with the observed stellar se- 
quence appearing in the BzK color region (see figure 44 
in the appendix). Thus, we can also use the BzK crite- 
rion of equation (1) as a robust method for star-galaxy 
separation. For sources detected with SN > 3, up to 
Kg = 20, we have excluded objects satisfying the follow- 
ing conditions; 1) equation (I) and 2) stcllarity parameter 
CLASS_STAR> 0.95 in SExtractor. In order to separate 
stars and galaxies for all of the ^'-detected objects up to an 
optical limiting magnitude z' = 26 that are fainter than the 
above Kg detected sources, we performed SED fittings for 
them with not only the BC03 models (Bruzual & Chariot 
2003) for galaxies but also the BPGS star atlas for stellar 
objects, including those corrected with Allen's Milky Way 
extinction curves. Wc extracted objects as stellar objects, 
which were fitted better by the latter than the former, 
from catalogue of z'-detected galaxies (see also subsec- 
tion 4.1). These stellar objects appear in the sequences of 
the BPGS star atlas in BBG color-color diagrams (see not 
only figure 44, but also figures 45, and 46 in appendix 1). 

4. Photometric redshifts 

4.I. Photometric and spectroscopic redshifts 

For all of the z'-dctcctcd galaxies remaining after the 
star-galaxy separation procedure as described in subsec- 
tion 3.4, we simultaneously obtained the photometric red- 
shift Zphot and the absolute magnitude Mn-opt correspond- 
ing to the stellar mass M, (see also appendix 7), by 
applying a photometric redshift code hyperz (Bolzonella 
et al. 2000). This is based on comparison of the ob- 
served ground-based u* BVRi' z' JKg photometry with a 
grid of stellar population synthesis models produced from 
Bruzual & Chariot (2003). From this BC03 model we 
prepared SED templates of a total of nine star-forming 
history models with exponential decaying star formation 
(SF) as SFR{t) oc exp(— t/r), with time scales for the star 
formation activity of r = 0.1, 0.3, 1,2, 3, 5, 15, and 30 Gyr, 
and we varied the age t from to 20 Gyr, and included 
the Constant Star Formation (CSF) model. All of the 
models are for solar metallicity. With the variations of 
E{B -V)= 0,0.1,0.2,0.3,0.4,0.5, and 0.6, we additionally 
applied Allen's Milky Way, Seaton's Milky Way, LMC, 
SMC, and Calzetti's extinction curves, to obtain the ex- 
tinction Ay. In selecting the best fitting one from five 
extinction curves for each object, firstly we determine a 
median photometric redshift Zm in five of them related to 




spec 

Fig. 1. Comparison between spectroscopic and photometric 
redsliifts for z'-detected galaxies with spectroscopic observa- 
tions. Filled diamonds and open circles represent objects with 
and without broad line emissions in their spectra, respectively. 
Dashed lines represent a deviation with Az/{1 + Zjpec) = 0.1. 

these extinction models, and secondly we have used one 
model with reduced minimization as the photometric 
redshift should be within 0.2(1 -fz^) around the median 
z„. We adopted a Salpeter IMF from 0.1 to 100 Mq 
throughout this work. 

In the AKARI NEP Deep field, wc also obtained spec- 
tra with the Keckll/DEIMOS, the Subaru/FOCAS, and 
the MMT/Hectospec for 420, 57, and 62 objects, which 
were selected from objects with R < 24, < 24, and < 
22 mag, respectively. The Kcckll/DEIMOS and the 
Subaru/FOCAS observations partially covered the NEP 
Deep field while the MMT/Hectospec observations cov- 
ered the NEP Deep field by two, 1 deg^ multi-object 
configurations. More details on the MMT/Hectospec ob- 
servations can be found in Ko et al. (2012). They are 
mostly detected as IR sources with the AKARI (Takagi 
et al. 2010). We have obtained 292 spectra for the z'- 
dctected galaxies, in which we can determine secure red- 
shifts for 230 objects are detected with more than two lines 
or clearly with [OH] line. They include 27 galaxies har- 
boring AGN with Broad Line (BL) emissions, which are 
spectroscopically classified as BL AGNs. Figure 1 shows 
the comparison between the spectroscopic and the pho- 
tometric redshifts of the 230 objects and their redshift 
discrepancy Az = Zphot — Zspec- The median accuracies 
of (Az) / {1 + Zspec) are 0.041 and 0.038 for the spectro- 
scopic samples without and with excluding the BL AGNs, 
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Fig. 2. Rest-frame color U — V for the spectroscopic sample 
and z'-detected galaxies. Top: Rost-frame Color Magnitude 
Diagram (CMD) of absolute V-band magnitude vs. U — V 
color. Objects in redshift ranges of 0.4 < z < 0.8, 0.8 < z < 1.2, 
and 1.2 < z are represented from the left to the right. Large 
symbols represent the spectroscopic sample. Red and blue 
symbols represent lower and higher SFR populations with 
sSFR < 0.1 Gyr-l and sSFR > 0.1 Gyr"! derived from 
the SED fitting, respectively. Magenta diamonds represent 
BL AGNs, in which diamonds in squares represent outlier BL 
AGNs with Az/(1 + Zspec) > 0.2 in the redshift comparison. 
Small gray dots represent the remainder in the z'-detected 
galaxies. Red and blue contours represent the distribution of 
lower and higher SFR population in the z'-detected galaxies 
with sSFR < 0.1 Gyr^i and sSFR > 0.1 Gyr^^, respectively. 
Dashed lines represent Bell's color boundary between star- 
forming and passive populations. Plain black crosses around 
U — V ^ —0.2 represent typical errors of estimated colors and 
mass. The vertical solid and steep slope dense dotted lines 
represent the limiting absolute V and (7-band magnitudes at 
the mean redshift, respectively. Open circles represent a pop- 
ulation categorized as type 1 QSOs at z = 0.2 — 0.7 in the 
catalogue of the SDSS/SWIRE sample by Hatziminaoglou et 
al. (2005). Bottom: Color distributions of the z'-detected 
galaxies. The color is scaled as {U - V) + 0.08(My + 20). 
Red and blue lines represent lower and higher SFR popula- 
tions with sSFR < 0.1 Gyr"! and sSFR> 0.1 Gyr-l derived 
from the SED fitting, respectively. 

respectively. In the spectroscopic sample, there are the 
22 and 52 outliers with Az/(1 + Zgpec) > 0.2 and > 0.1, 
in which 11 and 12 out of these outliers are BL AGNs, 
respectively. Thus, the photometric redshifts are consis- 
tent with their spectroscopic ones not only for the sam- 
ples except BL AGNs but also for even half of the BL 
AGNs. The outlier BL AGNs except one are confirmed as 
an Extremely Blue Object (EBO) in the rest-frame Color- 
Magnitude Diagrams (CMDs), which are selected with a 
criterion: 



{U-V)< -0.25(A/y + 22.0) + 0.7 , 



(2) 



as shown in figure 2. If EBOs are type 1 AGNs, this 
trend is reasonable since SED models of the AGNs were 
not included in SED fittings for the photometric redshift 
estimations. 

The AKARI MIR colors can alternatively extract can- 
didates associated with dusty AGN activities, in which 



we spectroscopically observed 14 agn-MbN3Rs and 2 s/a- 
MbN3Rs (see the details about them in subsection 6.1). 
We could confirm that 8 out of 14 agn-MbN3Rs and 2 
out of 2 s/a-MbN3Rs show the discrepancy with Az/ {1 + 
Zspec) > 0.1, respectively. All of these outliers are also 
identified as BL AGNs, in which 7 out of 8 the agn- 
MbN3Rs and 2 out of 2 the s/a-MBN3Rs are classified 
as the EBOs, respectively. As long as their redshifts are 
correctly obtained with spectroscopy, the overlap of the 
EBOs with the agn-MbN3Rs suggests that the dust emis- 
sion detected as the agn-MbN3Rs should coexist with op- 
tically blue emission classified as the EBOs in all of them. 
This supports not only the above expectations that AGNs 
are harbored in the EBOs and the agn-MbN3Rs but also 
a picture of a dusty torus surrounding AGN which can ex- 
plain the variation of their SEDs as frequently proposed 
in the AGN unified theory. 

We also checked the redshift discrepancy Az in various 
populations subclassified not only with the spectroscopy, 
but also with NIR and MIR photometry as N3Rs and 
MbN3Rs, which are summarized in appendix 3. Even 
though most of the spectroscopic samples are mainly se- 
lected from dusty populations with heavy extinction as 
MbN3Rs detected in the AKARI MIR photometry, their 
median accuracies of (Az)/(1 -I- z^pec) are less than 0.05 
for any subclassified populations of N3Rs and MbN3Rs 
except agn-MbN3Rs. Their photometric redshifts can be 
estimated well from the SED fitting with their major stel- 
lar emission of the host galaxies, in which the extinctions 
might not be a serious matter. Thus, the estimated photo- 
metric redshifts even for most of the galaxies can be used 
to reconstruct approximately their actual redshifts within 
Az/(1 + ) <0.05. 

4-2. Reconstructions of rest-frame magnitudes and col- 
ors 

From the derived photometric redshifts Zphot with the 
optical-NIR SED fittings, we can also estimate rest-frame 
magnitudes and colors of the z'-detected galaxies. For 
an object at a redshift z, a magnitude m\g(z) at a rest- 
frame wavelength Aq is approximately interpolated from 
the observed magnitudes as 

[A,, - A(z)] mi, + [A(z) - Afe] 



mxoiz) = 

A(z) = Ao(l + z), 



Xr — Xh 



(3) 
(4) 

where rrir and A^ (mf, and At,) are the observed magnitudes 
and central wavelengths of the band neighboring the red- 
side (blue-side) of the observing wavelength A(z). By us- 
ing the rest-frame magnitudes derived with Zphot, CMDs 
for the z'-detected galaxies can be obtained. The CMD 
fundamentally is almost the same as the color-mass dia- 
gram, which is a powerful tool for studying galaxy evolu- 
tion since it shows the bimodal galaxy distribution as early 
types concentrate on a tight red sequence while late types 
distribute in a blue dispersed cloud. This has been stud- 
ied not only in the local Universe with the SDSS (Blanton 
et al. 2003) but also in the distant Universe with surveys 
for z < 1 (BeU et al. 2004; Faber et al. 2007; Borch et al. 
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Table 3. The Parameters for IRBG Selections 
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2006) and for z > 1 (PanncUa et al. 2009; Brammer et al. 
2009). 

Figure 2 shows a CMD of the absolute l/-band mag- 
nitude Aly and U — V color. We display rest-frame 
^/-detected galaxies on the CMD, and excluded objects 
with less photometric accuracy of SN < 2 in the rest- 
frame V band. The rest-frame U — V color straddles 
the Balmer/4000 A break and is also essential to intro- 
duce classifications for their stellar populations in IRBGs, 
MbN3Rs, MmN3Rs, and BBGs as discussed in the follow- 
ing sections and the appendix. This interpolation scheme, 
applied for the z'-detected galaxies with z — Zphot, can be 
used to derive their My and rest-frame U — V colors. In 
figure 2, we can see that the distribution of galaxies is bi- 
modal up to z ~ 1 as previously shown (Bell et al. 2004; 
Faber et al. 2007; Borch et al. 2006). In figure 2, we can 
see that the red sequence and the blue cloud are separated 
with a criterion: 

{U-V) = 1.15-0.31z-0.08(Afy-51og/i-|-20) , (5) 

which is introduced for the definition of red sequence 
galaxies by Bell et al. (2004) and represented as dashed 
lines on the CMDs in figure 2. Thus, our catalogue of z'- 
detected galaxies with photometric redshifts is basically 
the same as those of other blank sky surveys. 

5. Classification with IR bump 

From the sample of z'-detected galaxies including the 
Balmer Break Galaxies (BBGs; see the details in ap- 
pendix 1), we extracted IR classified subpopulations with 
the AKARI/IRC photometry. The NIR IRC photometry 
can select IR bump Galaxies (IRBGs) which have a bump 
around 1.6 fim with the AKARI (Takagi et al. 2007) and 
the Spitzer (Berta et al. 2008). Both the IRBGs and the 
BBGs are characterized mainly by the emission from the 
stellar components. In this section, we will present three 
redshifted populations selected by criteria of the IR bump 
at 1.6 iim combining the Balmer break. In fact, the selec- 
tion with the 1.6 /xm can reduce the contamination from 
the low-z interlopers in the IRBGs at z > 0.4, which is also 
useful to classify MIR-bright populations. We also try to 
select AGN candidates as outliers with NIR colors redder 
than normal IRBGs. 

5.1. AKARI NIR color-color diagram and 1.6 fim bump 

A bump centered at 1.6 /im in the rest-frame is a major 
SED feature of galaxies except for the youngest ones for 
ages < 10 Myr, in which stellar components are dominated 
by type M stars. The 1.6 /zm IR bump is produced by the 
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Fig. 3. Two-color diagram of N2 - N3 vs. N3 - N4 in the 
field. Symbols of cross bar at (Af2 - iVS, Af3 - Ar4) = (0.1,0.5) 
represent typical errors of the color. Black line is superim- 
posed as a redshifted track of a BC03 model galaxy with 
T = IGyr and age= IGyr, on which the squares are related 
to the redshifts of z = 0.4,0.8, and 1.2 from lower-left to up- 
per-right. Other superimposed green, red, blue, cyan, and 
magenta lines are tracks of QSOl, TQSOl, BQSOl, QS02, 
and dusty torus IR templates, respectively. Redshift from 
z = to z = 2 corresponds to the upper-right to the lower-left 
along each line. Top: For objects detected in the N2,N3, 
and N4 bands. Arrows represent objects which are fainter 
than the threshold limit in the A'^2 or band. Red open 
circles and black crosses represent the objects at Zpf^at = 0.4 
and Zphot < 0.4), respectively. Bottom-Left: For the spectro- 
scopic samples. Blue, green, and red symbols represent the 
objects at 0.4 < Zspec < 0.8, 0.8 < Zspec < 1-2, and Zspec = 1-2, 
respectively. Diamonds represent BL AGNs. Bottom-Right: 
The same as on the left for the photometric samples. 



combination of the spectral peak in black-body radiation 
from low-mass cool stars as M stars and a minimum in 
the opacity of H" ion present in their stellar atmospheres 
with molecular absorptions (John 1988). Universally ap- 
pearing in most galaxies, the IR bump is one of the most 
useful SED features in their photometric measurements 
(Simpson & Eisenhardt 1999; Sawicki 2002). 

The AKARI NIR photometry in iV2iV3iV4 bands can 
trace the IR bump at z = 0.4 — 2. We have shown the 
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redshift tracks of the BC model galaxies with t = IGyr 
and age= IGyr on a two-color diagram of N2 — N3 vs. 
iV3 — in figures 3, where the solid squares on the tracks 
represent the redshifts at z = 0.4,0.8, and 1.2 from the 
bottom-left to the top-right. The tracks are mostly par- 
allel to the N2 — N3 axis up to z ~ 0.8, bending upward 
around z ~ 0.8, and are vertical along the N3 — NA axis at 
z ~ 0.8 — 1.2. This trend can be understood as the excess 
of redshifted IR bump is on the N2 side in the A^3 band at 
z ~ 0.4 - 0.8, the iV4 side in the iV3 band at z ~ 0.8 - 1.2, 
and in the NA band at z > 1.2. The redshift is the domi- 
nating factor in determining the AKARI NIR color tracks 
of galaxies without AGNs while the galaxy type is less 
effective. 

First, we extract A^3-detected ones as A^3 Red galaxies 
(N3Rs) possibly at z > 0.4 from the z'-detected galaxies 
with a criterion: 

N2~m> /jv23 ■ (6) 

Second, we can also subclassify the N3Rs into N23, N34, 
and N4 bumpers, which are defined as IRBGs of their IR 
bump detected with the A^2 and iV3 bands, the 7V3 and 
A^4 bands, and the A^4 band in three redshift ranges of 
z~ 0.4 -0.8, z~ 0.8 -1.2, and z> 1.2, respectively. The 
color criteria are represented as 

iV2 - 7V3 > /Ar23 n minAr34 < Ar3 - iV4 < maxAr34 , (7) 

where the parameters of /Ar23, niinAr34, and maxjv34 are 
given in table 3. 

The bottom-left of figure 3 represents a two-color dia- 
gram for the spectroscopic sample in the field, in which 
most of the galaxies at z^pec — 0.4 — 0.8,0.8 — 1.2, and > 1.2 
appear as N23, N34, and N4 bumpers, respectively. The 
bottom-right of figure 3 also represents another two-color 
diagram for the photometric sample in the field, in which 
the NIR color of galaxies at Zphot — 0.4 — 0.8,0.8 — 1.2, and 
> 1.2 distribute around the regions of N23, N34, and N4 
bumpers, respectively. 

When a galaxy harbors an AGN, however, IR emis- 
sion from the dusty torus around the AGN affect the NIR 
colors. Wc have also superimposed tracks of redshifted 
QSO SED templates in the SWIRE library three types 
of optically-selected type 1 QSOs (QSOl, TQSOl, and 
BQSOl) and two models of type 2 QSOs (QS02 and 
Torus). All of them show dominant emissions from the 
dusty torus at IR wavelength A > 1 /im. Thus, the NIR 
color criterion from equation (6) can be still effective for 
selecting not only galaxies at z > 0.4 but also AGN can- 
didates. Furthermore, all the QSO tracks are overlap- 
ping with the region of the N4 bumpers. In fact, most 
of the BL AGNs, selected with the spectroscopy, appear 
around the region of N4 bumpers. Thus, the N3Rs with 
A^3 — A^4 > —0.1 may include not only the N4 bumpers of 
IRBGs but also these AGN candidates, which will be also 
confirmed with the MIR classified AGNs in section 6. We 
win designate the N3Rs with iV3 - A^4 > -0.1 as N4 Red 
galaxies (N4Rs). 

* http://www.iasf-milano.inaf.it/ pollctta 
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5.2. N3Rs subclassified with Balmer break 




Fig. 4. The photometric redshift distributions for exclusively 
classified BBGs; e-BzKs, exclusively classified e-uRJs from 
the remainings after exclusion of the e-BzKs, and exclusively 
classified uVis from the remainings after exclusion of both 
the e-BzKs and the exclusively classified e-uRJs, from bot- 
tom to top. Shaded, single-hatched, and cross-hatched area 
represent p-, qp-, and s-BBGs. Left: For exclusively classi- 
fied BBGs and IRBGs. Dotted lines represent N23 bumpers, 
N34 bumpers, and N4Rs, from top to bottom. Right: For N3 
Red BBGs. Dashed, thin solid and thick solid lines represent 
the rest-frame V-detected galaxies, the N3Rs, all the N3 Red 
BBGs, respectively. 

As shown in figure 3, both spectroscopic and photo- 
metric redshifts of the N3Rs are almost consistent with 
the NIR color criterion 7V2 - iV3 > -0.3 of the IR bump 
at z > 0.4. As mentioned in subsection 4.1, the pho- 
tometric redshifts were derived from only the ground- 
based data without the AKARI NIR data, which means 
that the selection with the IR bump is independent of 
the photometric redshifts. Thus, the color criterion with 
A^2 — N3 > —0.3 is useful in order to exclude low redshift 
interlopers at z < 0.4 in the IRBGs as long as they are de- 
tected in the N2 (or K) and A^3 bands as discussed in sub- 
section 5.1. Furthermore, adding the AKARI/NIR pho- 
tometry to the ground-based optical- NIR photometry can 
salvage the BBGs in the IRBGs as shown in appendix 1.3. 
These are the merits in using the AKARI/NIR photome- 
try with the ground-based photometry. 

On the other hand, even though the IR bump criteria of 
N23 and N34 bumpers, and N4Rs correspond to three red- 
shift intervals as shown in figure 3, they cannot be robust 
for the redshift subclassification since the N4Rs possibly 
include AGNs as discussed in subsection 6.2. In redshift 
subclassifications, optical color criteria with the Balmer 
break arc more reliable than those of the A^3 — A^4 col- 
ors with the IR bump. As described in appendix 1, we 
have generalized the color criteria of BzKs at 1.4 < z < 2.5 
with BzK filter set (Daddi et al. 2004) to select BBGs at 
0.8 < z < 1.8 and 0.4 < z < 0.8 by using uRJ and uVi filter 
sets, called as uRJs and uVis, respectively. Furthermore, 
we have selected extended BBGs as e-BzKs and e-uRJs 
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Fig. 5. The same as figure 2 for N3 Red BBGs. The plots for 
0.4 < 2 < 0.8, 0.8<z< 1.2, and 1.2 < z from left to right. In 
general, colored open, solid, and circled symbols represent ob- 
jects classified as s-, p-, and qp-BBGs, respectively. Basically, 
red circles, green triangles, and blue stars represent objects se- 
lected with the two-color criteria as e-BzKs, o-uRJs, and uVis, 
respectively. Black-large open, solid, and circled symbols rep- 
resent mean values for the s-, p-, and qp-BBGs, respectively. 
Crosses on their black symbols represent their standard devi- 
ations. Only a part of samples arc plotted to reduce crowding 
in the diagram, whose fractions are presented as numbers in 
the right-bottom region of all figures. 

with the N3 band photometry which is the deepest one 
in the IRC bands (see appendix 1.3). However, these e- 
BzKs, e-uRJs, and uVis overlap each other. Thus, we 
introduced an exclusive subclassification of BBGs as first 
picked with the BzK criteria, second applied the uRJ cri- 
teria for those remaining after exclusion of the BzKs, and 
finally applied the uVi criteria for those remaining after 
exclusion of the uRJs and the BzKs. As shown on the left 
in figure 4, the redshift distributions of the exclusively 
classified BBGs except the uVis are similar to those of 
IRBGs as the uRJs and the BzKs correspond to the N34 
bumpers and the N4Rs, respectively. Most counterparts 
of the MIR bright populations at z > 0.4, detected with 
the AKARI MIR photometry in the field, are identified as 
the N3Rs. Thus, hereafter, we will take a subclassification 
of the N3Rs into the star- forming (s-), passively evolving 
(p-), and quasi passively evolving (qp-) BBGs with the 
exclusively subclassification. We will name these BBGs 
exclusively subclassified from the N3Rs as N3 Red BBGs 
(N3RBBGs); N3RBzKs, N3RuRJs, and N3RuVis. The 
redshift distributions of N3 Red BBGs are also shown on 
the right in figure 4. Figure 5 shows the CMD of My and 
[/ — y for these N3 Red BBGs. We can see that the sub- 
classification of the S-, p-, and qp-N3RBBGs is consistent 
with that using Bell's color boundary in figure 5. 

Since the iV3 band photometry is the deepest one in 
the IRC bands, as discussed in appendix 1.3, the N3 Red 
BBGs salvaged with the A^3 photometry are almost the 
same selected population as the rest-frame y-detected 
galaxies at z > 0.4. It can be confirmed as both of them 
show a similar redshift distribution as shown in figure 4. 
The rest-frame l^-detected galaxies are a kind of refer- 
enced sample as frequently plotted with the referenced 
contours on the CMDs in figures 2 and 5 (see also figures 
48 and 50 in the appendix). Thus, the N3 Red BBGs 
are suitable in the study for their stellar populations by 
using the CMD even in MIR bright populations as seen in 
section 6. 
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6. Classifications MIR SEDs 

The MIR multi-band photometry with the IRC has 
identified MIR-bright populations in the N3Rs as detected 
with more than an SN3 in MIR bands and classfied into 
starbursts or AGNs with two MIR color diagrams, which 
are called as MbN3Rs. The IRC MIR photometry can de- 
tect Poly-Aromatic Hydrocarbon (PAH) emissions from 
dusty star-forming regions for a subgroup in the IRC de- 
tected MbN3Rs. The subgroup with the PAH features can 
be classified with the MIR colors as starbursts, which are 
overlapped with the PAH luminous galaxies in some com- 
panion papers (Takagi et al. 2007, 2010). Another sub- 
group of the MbN3Rs without any PAH features or with 
weak ones is possibly AGN dominant. Some of them may 
overlap with Power-Law Galaxies (PLCs) and IR Excess 
Galaxies (IREGs), which are already classified with the 
Spitzer as AGN-dominant populations. The former has 
MIR SEDs weU fitted by a power law (Lacy et al. 2004; 
Stern et al. 2005; Alonso-Herrero et al. 2006; Donley et al. 
2007). The latter shows the large infrared-to-UV/optical 
flux ratio (Daddi et al. 2007; Dey et al. 2008; Fiore et al. 
2008; PoUetta et al. 2008), in which some may be iden- 
tified as sources with heavily obscured optical emissions. 
The MbN3Rs and their relatives at z > 0.4 detected in 
the field have typically TIR luminosities with > IO^^Lq, 
which correspond to those of LIRGs/ULIRGs as shown 
in appendix 4. In this section, thus, we will mainly study 
the MbN3Rs and their relatives by analyzing the MIR and 
NIR SEDs to reveal star forming and AGN activities of 
LIRGs/ULIRGs at z > 0.4 in the following sections. 

6.1. Starhurst and AGN in MIR bright NSRs 

The AKARI/IRC MIR photometry, in the S7, S9W, 
Sll, L15, L18Vt^, and L24 bands, is capable of tracing 
the features of PAH emissions at 3.3, 6.2, 7.7 and 11.3 ^m 
and Si absorption at 10 fiia (rest-frame) from dusty star- 
bursts in the MbN3Rs. The MIR two-color diagrams of 
S9W-SU vs. S'll-L15 {S9W/SU/L15) and S'll-L15 
vs. L15-L18W {Sn/L15/L18W) are useful to trace the 
redshifted PAH features as the strong PAH 7.7 fim emis- 
sion (the dip between PAH 6.2 and 7.7 /im) enters into 
the L15(S'll) band for higher-redshifted dusty starbursts 
in a redshift range of z = 0.4 — 2.6. In order to character- 
ize MIR SED features for the MbN3Rs, we present their 
observed colors on two-color diagrams of S9W/ S11/L15 
and S'11/L15/L18M^ as shown in figure 6. In order to 
classify them in the MIR two-color diagram, we have se- 
lected ~ 600 MbN3Rs from the N3Rs with the detection 
of S'iV > 3 in the MIR bands. 

Dusty starburst MbN3Rs (s-MbN3Rs) are distinguished 
from AGN dominated ones (agn-MbN3Rs) with solid lines 
in figure 6 as the boundary between the two populations, 
which are derived from the redshifted tracks of their MIR 
SEDs (sec figure 63 in the appendix). On the other hand, 
the agn-MbN3Rs arc selected with MIR color criteria in 
each of the redshift interval classes, which are 

(59W - 5*11) < 0.3 n (5*11- L15)> 0.0 (8) 
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Fig. 6. MIR color-color diagrams for the MbNSRs, which 
are also subclassified into s/p/qp-BBG/N3Rs and BL 
AGNs. Blue stars, green triangles, and red cir- 
cles represent s(p)-uVis, s(p)-(uRJs+uRN3s+BRN3s), and 
s(p)-(BzKs+uzN3s-|-BzN3s), respectively. Diamonds repre- 
sent the spectroscopic BL AGNs. Top-Left: The S9W - 511 
vs. 511 — L15 diagram for the observed MbN3Rs at 
0.4 < ^ < 0.8. Top-Right: The same as the Top-Left, for 
the MbN3Rs at 0.8 < ^ < 1.2 with the tracks in 0.8 < ^ < 1.4. 
Bottom: The 511 - L15 vs. L15 - LISW diagram for the 
MbN3Rs at 0.8 < ^ < 1.2 with the tracks in 1.0 < ^ < 2.8. 



at 2 = 0.4-1.0, 

{sm - sii) <{).?, n 

(S'll - Lib) < -1.5(5'11 - L15) + 1.2 
at 2 = 0.8-1.4, 

(S'll - L15) < 0.8 n (L15 - L18W) < 0.3 



(9) 



(10) 

at 1.0 < z < 2.8. Assuming that the MbN3Rs are at 
photometric redshift Zphot derived in subsection 4.1, the 
classification criteria at z = 0.4— 1.0, z = 0.8 — 1.4, and 
1.0 < z < 2.8 are apphed for the MbNSRs with the plroto- 
metric redshifts of 0.4 < Zp^ot < 0.8, 0.8 < Zphot < 1-2, and 
1.2 < Zphot, respectively, as shown in figure 6. The three 
redshift intervals roughly correspond to N3 Red uVis (N23 
bumpers), uRJs (N34 bumpers), and BzKs (N4Rs) in N3 
Red BBGs (IRBGs), respectively. Even with substitution 
of the photometric selections as the BBGs (IRBGs) for 
the redshift selections, the MIR two-color diagrams can 
be also applied to the distinction between dusty starbursts 
and AGNs. 

We have also fitted this AKARI/IRC photometry for 
the s-MbN3Rs observed at the 7 — 24 wavelength with 
the SED templates of starbursts (Siebcnmorgen & Kriigel 
2007) (hereafter S&K model, see also appendix 5.4). 
Even though we can roughly distinguish s-MbN3Rs and 
agn-MbN3Rs only with the AKARI MIR colors as their 
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Fig. 7. Rest-frame SED normalized at 1.6 fim for the 
MbN3Rs. Symbols are the same as figure 6 except for their 
colors. Blue, green, and red symbols represent the redshifts 
of 0.4 < 2 < 0.8, 0.8 < 2 < 1.2, and z > 1.2, respectively Dark 
and faint symbols represent AKARI IRC photometric data 
and ground-based photometric data, respectively. The av- 
eraged observed MIR SEDs are represented with solid lines 
in the figures on the right. Top-Left: For the sb-MbN3Rs 
in the spectroscopic sample. Superimposed lines represent 
dusty starburst SED templates of the S&K models used in 
the MIR SED fitting. Top-Right: The same as on the left for 
the sb-MbN3Rs in the photometric sample. Middle-Left: The 
same as the Top-Left for the s/a-MbN3Rs in the spectroscopic 
sample. Superimposed lines represent SEDs composed of the 
S&K models and the dusty torus model in the SWIRE SED li- 
brary with the ratio of 0.3 to 0.7 at 5 /im. Middle-Right: The 
same as on the left for the s/a MbN3Rs in the photometric 
sample. Bottom-Left: The same as the Top-Left for the ag- 
n-MbN3Rs in the spectroscopic sample. Superimposed lines 
represent the SWIRE SED templates of dusty torus, QS02, 
BQSOl, TQSOl, and QSOl from top to below at the long 
wavelength. Bottom-Right: The same as on the left for the 
agn-MbN3Rs in the photometric sample. 



MIR emissions are dominant from dusty starburst and 
AGN, some AGNs may still coexist with dusty starbursts 
as hidden in the s-MbN3Rs. In order to classify the 
s-MbN3Rs into two subgroups of starburst dominants 
and starburst/ AGN coexisting, we performed the MIR 
SED fitting with composite SED models, in which the 
dusty torus model in the SWIRE library is added to 
the dusty starburst S&K model with the mixing rates of 
10,20,30,40,50,60,70,80,90, and 100% in the rest-frame 
5-/um luminosity. We have selected sb-MbN3Rs and s/a- 
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Fig. 8. Correlation between the monochromatic luminosity 
vLy 7.7 at the rest-frame 7.7-/xm and the ratio of monochro- 
matic luminosities. Symbols are the same as figure 6 except 
for their colors and sizes. The blue, green, and red sym- 
bols represent the redshifts of 0.4 < 2 < 0.8, 0.8 < z < 1.2, 
and z > 1.2, respectively. Small symbols represent the ag- 
n-MbN3Rs classified with the IR SED fitting. Top; Ratio of 
the monochromatic luminosity vL^ lo at the rest-frame 10 /^m 
to the vL^ 7.7. Bottom: The ratio of the monochromatic lu- 
minosity vLi, 5 at the rest-frame 5 /^m to the vLu 7.7 for the 
MbNSRs. Left: For the sb-MbN3Rs. Right: For the s/a- and 
the agn-MbN3Rs. Large and small symbols represent the s/a- 
and agn-MbN3Rs, respectively. 

MbNSRs with the mixing rates of the dusty torus < 50% 
and > 50% as starburst-dominants and starburst/AGN 
coexisting, respectively, derived from the MIR SED fit- 
ting. The MIR SED fitting resuhs for the MbNSRs are 
consistent with their simple MIR color classification as 
MbNSRs fitted well with the high AGN mixing rate which 
overlaped with the agn-MbNSRs classified by the MIR 
color. Figure 7 shows the SEDs normalized at 1.6 /.tm 
for sb-, s/a-, and agn-MbNSRs where we had estimated 
the rest-frame monochromatic magnitudes (luminosities) 
with their spectroscopic and photometric redshifts. For 
the photometric sb-, s/a-, and agn-MbNSRs, their aver- 
aged MIR SEDs are also derived as shown with black lines 
on the right side. 

Rest-frame SEDs of the sb-MbNSRs, as shown in fig- 
ure 7, show deep dips around 5 and 10 /xm and a strong 
excess around 8 fiva. with weak ones around S and 11 /xm. 
The features agree well with SEDs of dusty starburst S&K 
models overplotted as thin lines. It indicates not only the 
consistency for this starburst classification with MIR col- 
ors and MIR SED fittings but also accuracy for their pho- 
tometric redshifts. On the other hand, the features are 
not so clear for the agn-MbNSRs. The 5 /im dip corre- 
sponds to the transition between the stellar emission and 



N2 - N3 redshift 

Fig. 9. Left: The same as figure 3 for the MbN3Rs. Symbols 
are the same as figure 6 except for their colors. Blue, 
green, and red symbols represent objects at the redshifts of 
0.4 < 2 < 0.8, 0.8 < 2 < 1.2, and z > 1.2, respectively. Right: 
Redshift vs. N3 - NA color. Top: For the sb-MbN3Rs. 
Middle: For the s/a-MbN3Rs. Bottom: For the agn-MbNSRs. 
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Fig. 10. The same as figure 7 for the MbN3N4Rs. Left: 
For the sb-MbN3N4Rs. Middle: For the s/a-MbN3N4Rs. 
Right:For the agn-MbN3N4Rs. 

the dust emission from their star-forming regions. The 
strong excess around 8 /im with weak ones around S and 
11 /im and the dip around 10 /im can be also naturally 
explained with the PAH 7.7, S.S, and 11. S /im emissions 
and Si 10 /im absorption in typical MIR emission from 
their dusty star-forming regions. 

Figure 8 shows similar trends to those which appeared 
in figure 7 as the dips at rest-frame 5 and 10 /im are deeper 
in the sb- MbNSRs than in the s/a- and agn-MbNSRs. The 
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trends can be also naturally explained as MIR SEDs of 
s/a-MbN3R consist of a smooth power-law or convex com- 
ponent adding to an original component with deep con- 
cave dips at 5 and 10 fim^ which agree well with the super- 
imposed SEDs composed of dusty starburst S&K models 
and the dusty torus model in the SWIRE SED library 
with the ratio of 0.3 to 0.7 at 5 /xm. The composition of 
starburst and dusty torus is typically characterized by the 
monochromatic luminosity ratio of s/a-MbN3Rs shown in 
figure 8. This is consistent with a picture of an obscured 
AGN coexisiting with dusty star- forming region in an s/a- 
MbN3R which contributes to the additional smooth com- 
ponent and the emission with the 5 and 10 fim dips in the 
MIR SED, respectively. 

6.2. NIR colors of MbNSRs 

The NIR color diagram of figure 9 for the MbN3Rs can 
also support the existence of the dusty torus hidden in the 
agn- and s/a-MbN3Rs with the excess in the iV4 band. 
The agn- and s/a-MbN3Rs even at z < 1 appear in the 
region of A^3 — > —0.1, which is the same trend as 
for the BL AGNs and expected from the tracks of the 
dusty torus and QSO models as shown in figure 9. For 
stellar compornent with the IR bump, the mean iV3 — NA 
color track, represented as solid lines in the right plots 
of figure 9, can be roughly approxinated as a function 
N3NA{z) of z: 



N3NA{z) 



0.74Z-0.84 (0.4<z<1.2) 
0.048 (z>1.2) 



(11) 



Thus, we have selected AGN candidates with an NIR color 
criterion: 



7V3 - iV4 > N3NA{z) + AiV34 , 



(12) 



where we have taken a typical iV3 — A^4 color error A{N3 — 
iV4) into account and introduced AN 34 = 3A(iV3 - NA). 
The galaxies, selected with this additional NIR color cri- 
terion from the N3Rs, will be called N3N4 Red galax- 
ies (N3N4Rs). The boundary of the NIR color criterion 
is shown as dotted dash lines in the right plots of fig- 
ure 9, where A{N3 - iV4) = 0.1 as a typical N3 - NA 
color error for the MbN3Rs, which selects MIR bright 
N3N4RS (MbN3N4Rs). As shown in figure 10, the SEDs of 
agn- and s/a-MbN3N4Rs are similar to spectroscopically 
confirmed BL AGNs in agn-MbN3Rs. The redder color 
N3 - Ni in the agn- and s/a-MbN3N4Rs is consistent 
with the existence of a dusty torus, which dominates stel- 
lar components of the host. On the other hand, the SEDs 
of sb-MbN3N4Rs are still similar to those of dusty star- 
bursts as most of the sb-MbN3Rs. The sb-MbN3N4Rs, 
even with the redder color A^3 — A^4, may be still dusty 
starburst dominant. Thus, the MIR SED classification for 
the sb-MbN3Rs is robust to distinguish dusty starbursts 
from AGNs while the NIR color classification can be ef- 
fective only to select bright obscured AGNs. 

6.3. MIR marginally- detected NSRs 

We should remark that the MIR multi-band photome- 
try with the IRC can select unique samples for studying 
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Fig. 11. The same as figure 6 for the 
pP-MmN3Rs. Dark and faint symbols represent 

the sb- and non sb- pP-MmN3Rs, respectively. 
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Fig. 12. The same as figure 11 for the fP-MmN3Rs. Dark 
and faint symbols represent sb- and non sb-fP-MmN3Rs, re- 
spectively. 

star formation and AGN activity as not only the MbN3Rs 
but also sources detected only in some of the S7, S9W, 
Sll, L15, L18W, and L24 MIR bands. As shown above, 
the Sll, L15, L18W, and L24 of the IRC are effective in 
detecting the 7.7 /xm PAH emission from dusty starbursts 
at z > 0.4. Even though we have selected 600 MbN3Rs 
with the detection of SN > 3 in more than two MIR bands 
in subsection 6.1, we have also excluded MIR sources 
in the N3Rs as detected only in one or two MIR bands 
with detection of SN > 3 and in other MIR bands with 
detection of SN < 3, which are called MIR marginally- 
detected N3Rs (MmN3Rs), whose MIR magnitudes and 
colors can be still estimated as the MmN3Rs are mimics of 
the MbN3Rs. We selected --600 MmNSRs, which are sub- 
classified into ~320 possible PAH emitting (pP-)MmN3Rs 
and --270 faint PAH emitting (fP-)MmN3Rs. The former 
can be PAH 7.7 /im emitters as detected in the Sll, L15, 
and L18 bands at z = 0.4 — 0.65,0.65 — 1.5, and > 1.5, re- 
spectively, while the latter even have MIR emission but 
possibly not related to the PAH emissions. Thus, the pP- 
MmN3Rs and fP-MmN3Rs could be candidates for dusty 
starbursts and AGN mixtures or less-active star-forming 
with and without the PAH emissions in their selection 
criteria, respectively. 

Following the MIR color classifications in subsec- 
tion 6.1, we checked the difference between the pP- and 
fP-MmN3Rs on the MIR color-color diagrams as shown 
in figures 11 and 12. Even though the errors in the MIR 
colors of the MmN3Rs are large, most of the pP-MmN3Rs 
appear in the area of dusty starburst in the MIR color- 
color diagrams, which are similar to the s-MbN3Rs. On 
the other hand, the fP-MmN3Rs are more dispersed on 
the diagrams, which suggests that the fP-MmN3Rs mainly 
consist of star-forming and AGN mixtures. In order to 
classify the MmN3Rs into two subgroups of starburst and 
non starburst, we also performed the MIR SED fitting 
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Fig. 13. The same as figure 7 for the MmNSRs. Left: For 
the sb-MmN3Rs. Right: For the non sb-MmN3Rs. 
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Fig. 14. The same as figure 8 for the sb- and non 
sb-MmN3Rs. 

with the same composite SED models as described in sub- 
section 6.1. We selected sb- and non sb-MmN3Rs with 
the mixing rates of the dusty torus < 50% and > 50% 
in MIR SED fittings as starburst-dominants and possi- 
bly non starburst, which are represented as dark and thin 
symbols, respectively, in figures 11 and 12. Then, we cate- 
gorized these starburst dominant MmNSRs (sb-MmN3Rs) 
as mimics of sb-MbN3Rs. 

In fact, - 240 out of - 320 pP-MmN3Rs are sb-pP- 
MmN3Rs, while only - 110 out of - 270 fP-MmN3Rs 
are sb-fP-MmN3Rs. Thus, the classifications from the 
MIR colors for pP-MmN3Rs and fP-MmN3Rs arc roughly 
consistent with those from the MIR SED fittings for sb- 
MmN3Rs and non sb-MniN3Rs. Figure 13 also shows av- 
eraged MIR SEDs of sb-MmN3Rs and non sb-MmN3Rs. 
We can see that the averaged MIR SED of the sb- 
MmN3Rs have dusty starburst features with deep dips 
around 5 and 10 /im and excesses around 8 and 11 fim in 
the SEDs, which arc similar to the sb-MbN3Rs as shown 
in figure 7. We can also sec these features in the lumi- 
nosity ratios of vL^^io/i/L^j^r and vL^^^/vL^^t/j as shown 



Fig. 16. The same as figure 7 for the MmN3N4Rs. 

in figure 14. On the other hand, the averaged MIR SED 
and the luminosity ratio of the nonsb-MmN3Rs are sim- 
ilar to those of the s/a-MbN3Rs. We also tried to find 
AGN candidates hidden in the MmN3Rs with the NIR 
color criterion of equation (12) as applied in the MbN3Rs 
in subsection 6.1. With A(A^3 — A^4) = 0.15 as a typical 
7V3 — NA color error for the MmNSRs, the boundary of 
the criterion is represented again as dotted-dashed lines 
in figure 9, which select MmN3N4Rs. As shown in fig- 
ure 16, the sb- and nonsb-MmN3N4Rs are similar to the 
sb- and agn-MbN3N4Rs, respectively. 

Thus, the distinction between starbursters and AGNs 
with features in MIR/NIR SED is basically valid even 
in the MmNSRs the same as the MbNSRs. Thus, we can 
treat sb-MmNSRs as mimics of s- MbNSRs in the following 
sections. 

6.4. MIR faint NSRs and BBGs 

From the opposite point of view to subsections 6.1 
and 6.3, weak emission or non-detection in these MIR IRC 
bands are still important constraints for selecting galax- 
ies with low extinction or AGNs with weak dusty torus 
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emission. These MIR faint sources, having less than the 
hmiting magnitudes of those in table 2 in all the MIR 
bands (S9, Sll, L15, L18W, and L24), can be classi- 
fied into a subclass population of N3Rs, which we call 
as MIR faint N3Rs (MfNSRs) as alternatives to not only 
the MbNSRs but also the MmN3Rs. The MfNSRs may 
in general also consist of star-forming and passive popula- 
tions in the BBGs with weak dust emission. For example, 
passive uVis, uRJs, and BzKs without the 7.7 /im PAH 
emission at 0.4 < z < 0.8, 0.8 < z < 1.2, and 1.2 < z are 
expected to be faint in the Sll, S11/L15, and L15/L18W 
bands, respectively. Thus, the MfNSRs can be subclassi- 
fied into s-and p-BBGs with the two-color diagram. 

The MfNSRs at z > 1.2 form a limited sample in our 
shallow survey depth in the J, K and A^S bands, and 
they are mostly identical to MIR faint BBGs (MfBBGs), 
MfNSRs at z < 1.2 are a part of MIR faint uVis (MfuVis). 
In order to study the evolution and mass dependence in 
star formation, these MfNSRs and MfBBGs are also im- 
portant as reference samples compared with the MbNSRs 
and MmNSRs for analyzing their evolutionary features of 
SFR, extinction, and in section 8. We will sometimes 
refer to MfBBGs as MfBBGs not classified as MfNSRs in 
the following. 

6.5. Comparison with X-ray and Radio Observations 

Before closing this section, we will quickly compare 
the MIR SED classification with recent results at X- 
ray and radio wavelengths. AGNs being harbored in 
the agn-MbNSRs and the s/a-MbNSRs are also indicated 
from a preliminary look at our recent Chandra X-ray 
Observatory (CXO) data on this field (Krumpe et al. in 
prep.; Miyaji et al. in prep), which show that ^^50% of the 
agn-MbNSRs and -15% of the s/a-MbNSRs have coun- 
terparts of X-ray sources. Thus, the trend for the X-ray 
sources also support that the MIR SED classification is ef- 
fective for selecting AGNs. Furthermore, as shown in fig- 
ure IS, the nonsb-MmNSRs are possibly harboring AGNs 
with dusty starbursts. If it is true, their X-ray emission 
can be detected with stacking analysis of their CXO data. 
Thus, the CXO data can be effective to study AGN activ- 
ities not only for the MbNSRs but also for the MmNSRs. 

We have also made a radio observation with the WSRT 
at 1.4 GHz of ~ 1.7 deg^ covering the field, which found 
~ 500 sources with > 0.1 mJy in the field (White et al. 
2010). The radio source counts are consistent with a 
two population model with radio loud sources and ra- 
dio faint sources with < 1 mJy, in which it was assumed 
that the former and the latter are powered by AGNs and 
star-forming galaxies as in previous studies of other deep 
fields. The radio loud and the sub-mJy populations pos- 
sibly overlap with the agn- and s-MbNSRs, respectively. 
However, as the s-MbNSRs were subclassified to the star- 
burst dominates (sb-MbNSRs) and the starburst/AGN 
co-cxistence(s/a-MbNSRs) in subsection 6.1, the sub-mJy 
radio populations can be also subclassified to starburst 
dominant and starburst/AGN co-existence populations. 
Thus, it will be interesting to compare the MIR popu- 
lations with the X-ray and radio sources, which will be 



studied in a future paper. 

7. Calorimetric studies vi^ith MIR SEDs 

For the s-MbNSRs and the sb-MmNSRs classified as 
dusty starbursts, their MIR SEDs with PAH emissions 
and Si absorption should contain rich information not only 
as distinguished from AGNs but also about the physical 
properties in their star-forming regions such as their bolo- 
metric luminosities, SFRs, and extinction, which were de- 
rived by calorimetric schemes in MIR SED analysis. 

7.1. Bolometric Luminosity and SFR 

As long as AGN activity is negligible in a galaxy, young 
stars are dominant sources for the emission. Thus, the 
contribution of SFRs are nearly equal to the bolomet- 
ric luminosity for star-forming galaxies. We can estimate 
the total emitting power with an SED analysis for multi- 
wavclcngth photometric data and convert it to the SFR. 
SFRs are derived not only from intrinsic UV continuum 
luminosity by assuming an extinction law for all the z'- 
detected galaxies, but also from IR luminosity for the s- 
MbNSRs and the sb-MmNSRs. 

When massive stars are not obscured by dust, SFR on 
a star-forming region can be directly traced with their 
UV continuum luminosity in the rest-frame wavelength 
ISOOA. Thus, the UV luminosity can be translated into 
the SFR by the relation: 



SFRuv [Meyr-i] 



i^(1500A) 



8.85 X 1027[erg s-^Hz" 



(13) 



When massive stars are born in dust obscured star- 
forming regions, their UV radiation is reprocessed into 
IR emission. In this case, the SFR can be traced mainly 
with their TIR luminosity Lm. Even without using their 
TIR luminosity, as frequently applied, the intrinsic UV 
continuum luminosity Luy-cor can be reconstructed with 
a correction for the classical dust extinction Auv-.opt from 
the observed UV luminosity Luy-obs (sec also the details 
in subsection 7.S), and derive SFRuv ^cor from Ljjv.cor 
with equation (IS) and/or directly optical- NIR SED fit- 
tings with synthesis models. The scheme for SFRijv;cor 
can be applied even for all the z'-detected galaxies includ- 
ing MIR faint galaxies. 

For the s-MbNSRs and the sb-MmNSRs, on the other 
hand, with integrating a fitted model SED from 8 to 1000 
^m in the rest-frame with correcting for the photometric 
redshift Zphot , we can obtain the total IR (TIR) luminosity 
Ljft of the emission from the dusty star forming regions. 
According to Kennicutt (1998), thus, the SFR can be es- 
timated from Lifi as 

SFRMy^syr-^] = ,_,,,,%^,-.y (14) 

where we assumed the Salpeter IMF in the optical-NIR 
SED fitting as shown in subsection 4.1. In general, the 
total SFR in a galaxy is derived as a sum: 



SFR 



IR+UV 



SFR 



IR 



SFR, 



UV;obs , 



(15) 
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since dust can partially obscure galaxies. It is noted that 
SFR;7y;obs is derived from the observed UV luminosity 
LuV;obs without extinction correction. 




i'L^[7.7^m)/L0 iyL^(7.7^m)/LQ 




^0° io'° lo" -lo' io'° lo" 

v\ (7,7nm)g,/Lg v\ (7.7/im)sj/L3 



Fig. 17. Rest-frame 7.7-/im monochromatic luminosity 
t^I^iy 7.7 ij.m vs. the TIR luminosity Lju derived from the IR 
SED fittings for the s-MbN3Rs. Symbols are the same as in 
figure 6 except for their colors and sizes. Blue, green, and red 
symbols represent objects at the redshifts of 0.4 < z < 0.8, 
0.8 < z < 1.2, and z > 1.2, respectively. Large, medium, and 
small symbols indicate the results which fit well with the 
S&K model SEDs of Ay-sK = 2.2,6.7, and 17.9, respectively. 
Top: The L/ij derived without excluding the contribution 
from the hidden AGN. Bottom; The Lm derived with 
excluding the contribution from the hidden AGN. Left: For 
the sb-MbN3Rs. Right: For the s/a-MbN3Rs. 

7.2. Total IR Luminosities from MIR SED 

By means of the S&K models as a library of SED tem- 
plates required in MIR SED fittings, we could derive TIR 
luminosity from the fitted IR SEDs for the s-MbN3Rs. 
Thus, we could obtain a relation between the rest-frame 
monochromatic 7.7 fim luminosity vL^ti and the TIR 
luminosity Lm for all the s-MbN3Rs including both sb- 
MbNSRs and s/a-MbN3Rs as shown in figure 17. The 
top and bottom diagrams in figure 17 indicates the results 
in the Lm. estimation without and with the subtraction 
of the obscured AGN contribution, respectively. Before 
the subtraction of the AGN contribution, TIR luminos- 
ity range of the s/a-MbN3Rs is larger than those of the 
sb-MbN3Rs as shown in the top diagrams in figure 17. 
After the subtraction, the discrepancy in TIR luminosi- 
ties between the sb- and s/a-MbN3Rs was reduced. This 
suggests that the subtraction is efi'ective to derive TIR lu- 
minosities even for the s/a-MbN3Rs, which is also essen- 
tial to estimate their SFRs as discussed in subsection 8.2. 

Even though the TIR luminosity correlates well with 
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Fig. 18. Relation between TIR luminosity L/^ estimated 
from the S&K model and the rest-frame luminosities of 
Lv;7.7 fim and L^;io (jm- Symbols are the same as in fig- 
ure 6 except for their colors and sizes. Blue, green, and red 
symbols represent objects at the redshifts of 0.4 < z < 0.8, 
0.8 < z < 1.2, and z > 1.2, respectively. Large, medium, and 
small symbols indicate the results which fit well with the S&K 
model SEDs of Ay-SK = 2.2,6.7, and 17.9, respectively The 
line represents the relation approximated with equation (16). 

the rest-frame 7.7-fim monochromatic luminosity vL^i ^ 
as shown in figure 17, this estimation of the TIR lumi- 
nosity with the MIR SED fitting seems to be still a kind 
of extrapolation to obtain the Far InfraRed (FIR) emis- 
sion around ~ 1000 /ini from the NIR/MIR side only up 
to ~ 24 /.tm with the best-fitted model SED. However, we 
would like to remark that it can be more than an extrap- 
olation scheme since the AKARI MIR multi-band pho- 
tometry can certainly trace not only the PAH emission 
but also the Si absorption feature as already shown in 
figures 7 and 8. The absorption feature depends on the 
optical depth in galaxies and corresponds to the TIR lu- 
minosity as discussed in the following. 

We can see multi-sequences in figure 17. The splitting 
into the multi-sequences is mainly caused by the differ- 
ence in the "mean" optical depth, which can be charac- 
terized as an SED model parameter of "visual extinction" 
Av-SK at the edge of dusty starburst nucleus in the S&K 
model, which is indicated by the symbol size. In MIR 
SED fittings with the S&K model library, we have taken 
the model SEDs with Ay^sK = 2.2,6.7, and 17.9 corre- 
sponding to large, medium, and small symbols in figures, 
respectively. Even though Ay-SK is not equivalent to the 
observed visual extinction Ay assuming spherical geome- 
try for star forming regions in the S&K model, this param- 
eter Av-^sK is quantitatively related to the optical depth 
of the observed dusty starbursts as the model MIR SEDs 
with larger Ay-^SK show a deeper Si 10 /im absorption 
feature since it is induced with the dust self-absorption 
process in dense dusty regions. In fact, the s-MbN3Rs, 
with not only PAH emissions but also Si 10 /im absorp- 
tion, have various ratios of vL^.^iq/vL^.j ^ between 10 and 
7.7 ^la. as already shown in figure 8. Then, the observed 
ratio vLu-^iq/vL^-tj contains the information about the 
"mean" optical depth related to the denseness and the 
size of a dusty star-forming region, which is parametrized 
with the Ay-sK in the S&K models (see the details in 
appendix 5.4). 



No. ] Dusty IRBGs/LIRGs at z = 0.4 - 2 

Table 4. Fitting results with equation (21). 
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(RMS) 
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3.3 
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0.2 


1 SI 




fixed 1 


1.4 
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0.1 


1.69 


s/a-MbN3R 


0.41 ±0.08 


3.2 


± 


0.3 


1.47 




fixed 1 


1.0 


± 


0.1 


1.74 


sb-MmN3R 


0.22 ±0.09 


4.6 


± 


0.3 


2.34 




fixed 1 


2.0 


± 


0.1 


2.57 


non sb-MmNSR 


0.29 ±0.12 


3.7 


± 


0.4 


2.46 




fixed 1 


1.5 


± 


0.1 


2.68 



Values are means and ± standard deviations. 



Tfius, the MIR SED fitting with the S&K model can 
analyze a variational degree of "mean" optical depth pos- 
sibly related to the Si absorption depth. The "mean" 
optical depth is another infiuential parameter for estimat- 
ing the TIR luminosity from the MIR SED fitting. Thus, 
it is expected that their TIR luminosity is correlated not 
only with the rest-frame monochromatic 7.7 /im PAH lu- 
minosity vhy^-j 'j but also with the Si absorption depth at 
least for the selected starburst samples of the s-MbN3Rs 
and the sb-MmN3Rs as 



i^IR;SK 10 X 



vL 



X vL 



1/7.7 



(16) 



where the units of Lif(,vLi,-; -! and vLi/iq is [erg s~^] and 
a ~ — 1. This approximate formulation can work well as 
shown in figure 18. 

Even though the S&K models can reproduce MIR SEDs 
of the s-MbN3Rs and the sb-MmN3Rs well, they are pro- 
duced with some simplifying assumptions such as a spher- 
ical geometry with constant dust distribution for the mod- 
eled starburst region. Even though we could confirm the 
accuracy of the estimation for the TIR luminosity with 
complementary observations at longer wavelengths than 
MIR as FIR and submillimeter as remarked in section 1 1 , 
we have not yet obtained them in the field. At this mo- 
ment, we can not rule out a possibility that the TIR lu- 
minosity Lifi-sK, estimated with the fitted S&K models, 
includes a systematic offset. Thus, we will introduce a 
factor fsK in the conversion from the model estimated 
LiR-sK to the real Lir: 



Lib, — /skLih-sk 



(17) 



As long as the proportional factor fsK — 0.3 — 1, the es- 
timation of Ljii is consistent with previous work. For 
example, taking fsK — 1, we can roughly reproduce a re- 
lation Ljfi — 17.2L5_8.5 ^ni by Elbaz et al. (2002). We 
found an offset between the S&K model and the C&E 
model, in which the TIR luminosity with the C&E model 
can be reproduced with fsK — 0.3 (see figure 64 in the ap- 
pendix). In the following, we take fsK = 1 as the default 
value. 
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Fig. 19. Correlation between two extinctions at 1500A de- 
rived from the optieal-NIR SED fitting and the luminos- 
ity ratio of the total to observed UV as Ai/y-apt vs. the 
SFRijij^UV / SFRuY^obs- Symbols are the same as in fig- 
ure 18. Dashed line represents a proportional law with equa- 
tion (21), obtained with fitting the sample. Left: For the 
sb-MbN3Rs. Right: For the s/a-MbN3Rs. 




Fig. 20. The same as figure 19 for the MmNSRs. Dashed line 
represents a proportional law with equation (21), obtained 
with fitting the sample. Left: For the sb-MmN3Rs. Right: 
For the nonsb-MmNSRs. 

7.3. Classical and calorimetric extinctions 

In this subsection, we will derive extinctions at UV 
wavelength ~ 1500A with two schemes and compare them. 
One scheme is with optical-NIR SED analysis, deriving 
classical UV extinction Aijv,opt related to common visual 
extinction Ay. Another is with the IR SED analysis, de- 
riving calorimetric extinction Ajjv.iR as a ratio of total 
emission to observed UV emission. 

For all BBGs/IRBGs including galaxies even without 
MIR detections, visual extinction Ay can be derived from 
the multicolor optical-NIR SED fitting with stellar pop- 
ulation synthesis models as the BC03 model. In a dust- 
screen model, the dust extinction magnitude A\ is pro- 
portional to the optical depth at a wavelength A to char- 
acterize the difference between intrinsic flux i^i„t(A) and 
observed flux Foh<,(A) as 



,(A) = i^™t(A)10- 



(18) 



With selecting a reddening curve of fc(A) from extinction 
models: the Milky way, the LMC, the SMC and Calzetti 
law, extinction A\ at a wavelength A can be represented 
with the visual extinction Ay or the color excess E[B — V) 
as 
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Fig. 21. Extinction difference AAijy vs. the 
monochromatic luminosity ratio of vLi,n^ /vL^';,-j , 
where /^Auv = Auv;lR - Auv;opt is de- 
rived from the calorimetric estimated extinction 
AuV;lR = 2.5 \og{SFRiii+uv/SFRuv:obs) and the 
classical extinction Any-apt ■ Symbols are the same as in 
figure 18. 



A^ 



fc(A) 



Av ^k{X)E{B -V) . 



(19) 



where Ry = 3.1,2.72, and 4.05 is for both the Milky way 
and the LMC, for the SMC, and for the Calzetti law, re- 
spectively, and the reddening curves are normalized as 
k{B) — k{V) = 1. In the following, we have taken the 
best-fitted reddening curve in the optical-NIR SED fit- 
tings for each galaxy to convert the obtained Ay^opt to 
Auv;opt = ^1500 at A =1500A as a typical "optically ob- 
served" UV extinction. 

Even though the above estimation of Av;opt has been 
frequently applied in studying extinction of distant galax- 
ies, it should be recalled that it is based on a simple 
dust-scrcen model, which cannot guarantee the correct 
description of the radiation fields in nuclear starbursts 
expected in the s-MbN3Rs and the sb-MmN3Rs. Thus, 
we also tried to estimate calorimetric extinction AjjvjR 
for the s-MbN3Rs and the sb-MmN3Rs. If the dust and 
stellar components are mixed well in a galaxy as nu- 
clear starbursts in local LIRGs/ULIRGs, extinction es- 
timations with the IR SED analysis may be more ac- 
curate than the common scheme for the classical extinc- 
tion Auv;opt from visual extinction Ay. Total bolometric 
luminosity Lm^uy from the star-forming regions in s- 
MbN3Rs and sb-MmN3Rs can be estimated as Lju^jjy — 
Liu + Luy-obs-, which is the sum of the TIR luminosity 
L/zj derived from the MIR photometric SED analysis in 
section 7 and the observed UV luminosity Luy-obs not cor- 
rected for dust extinction derived from optical-NIR pho- 
tometric SED fitting. Using Ljfj^jjv and Lijy-obs, we can 
derive an alternative extinction with a calorimetric scheme 



A 



UV-IR ' 



2. 5 log 



■•IR+UV 



UV;obs 



(20) 



where we have taken the SFR ratio as the luminos- 
ity ratio {L I R+uv / Luy-obs) — {SFRm+uy / SFRuy-obs)- 
Figures 19 and 20 show the comparison between the 
calorimetric extinction Auy-jfj and the classical extinc- 
tion Ajjy-^opt = ^1500 for the s-MbN3Rs and the sb- 



MmNSRs, respectively, where we have taken fsK = 1 in 
the estimation of Lju. We can see that the calorimet- 
ric extinction Ajjy-jji with the IR SED analysis shows 
a proportional trend for the classical extinction Ajjy-opt, 
which suggests both extinctions are different, yet almost 
the same. 

However, there are some offsets between Ajjy-jft. and 
Auv ■opt as shown in figures 19 and 20. We have fitted 
the trend with a form as 



hjVJR — O-AAuy-^opt + AAf/y 



(21) 



As the fitting results arc summarized in table 4, mean off- 
sets of sb-MbN3Rs, s/a-MbN3Rs, and sb-MmN3Rs are 
significant compared with the RMS of residuals while 
those of non sb-MmN3Rs are not. The result suggests that 
the offsets of sb-MbN3Rs, s/a-MbN3Rs, and sb-MmN3Rs 
may be systematic ones. With fixing a^i = 1 as a propor- 
tionality between Ajjy-jR and Ajjy-^opt assumed for sb- 
MbN3Rs, s/a-MbN3Rs, and sb-MmN3Rs, the systematic 
offsets are (AAjjy) = 1.4,1.0, and 2.0, respectively, which 
are represented with a dashed line as shown in figures 19 
and 20. There are two ways to understand the system- 
atic offset; 1) the Auy-m with the S&K model tends to 
overestimate, or 2) the Ajjy-^opt with the classical extinc- 
tion laws tends to underestimate. At least for s-MbN3Rs, 
the former recommends to take fsK — 0.3 while the latter 
suggests fsK — 1- On-going or near-future observations 
in submillimeter and FIR wavelengths can determine the 
parameter fsK as discussed in section 11. 

The MIR SED fitting with the S&K model derives 
not only Ajjy-jR from L^-^sK but also visual extinction 
Ay-sK in the model parameters. In figures 19 and 20, 
the smaller and larger symbols represent galaxies fitted 
with larger and smaller Ay-sK, respectively. We can see 
another trend as sb-MbN3Rs with a larger Ay-sK has a 
larger offset in extinction /S.Aijy due to the variation of 
Auv;iR even in the same Ajjy-opt sample. This varia- 
tion of Ajjy-jR, can be caused by that of interstellar den- 
sity, which also determines a parameter Ay-sK in the 
S&K models. As discussed in subsection 7.2, the Ay-^sK 
characterizes the monochromatic luminosity ratio between 
rest-frame 10 and 7.7 fmi;i'L^io/i'L^Y-7 as the effect on 
"mean" optical depth in Lm/vLj^-j -j represented in equa- 
tion (16) (see also figures 17 and 18). Then, it is ex- 
pected that we can see some correlations between /S-Ajjy 
and vLiJio/uLijjj. In figure 21, indeed, we can see that 
sb-MbN3Rs with a larger Ay-sK have not only a larger 
/S-Ajjy-jR, but also a deeper 10-/im decrement, which may 
be caused by the Si self-absorptions as remarked in sub- 
section 6.1. 

These trends can be naturally explained by different 
characteristics of the extinctions of Ajjy.opt and Auy-m 
as follows. The calorimetric extinction Auy-jji may trace 
the dust extinction in dense dust-obscured areas in star- 
forming regions as the dust column density becomes large 
enough to cause the Si self-absorption. In contrast, the 
classical extinction Ajjy-^opt tends to trace relatively lower 
optical depth areas such as optical emissions which are 
mainly detected from the lower optical depth areas. Thus, 



No. 



Dusty IRBGs/LIRGs at z = 0.4 - 2 
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the extinction difference AAjjv may be related to the ge- 
ometrical variations in their star-forming regions and dust 
distributions in the s-MbN3Rs and sb-MmN3Rs. 

8. Evolution of SFRs and extinctions 

As examined in section 7, we could derive the TIR lumi- 
nosity and the extinction for s-MbN3Rs and sb-MmN3Rs 
at z = 0.4 — 2, which can be applied to estimate the SFR. In 
this section, we will reconstruct their evolutionary features 
and mass dependence of extinction and SFR in galaxies 
up to z = 2. 




Fig. 23. The same as figure 
For Av;iR and SFRir+uv- 

SFRuY-^arr- 



22 for sb-MmN3Rs. Left: 
Right: For Av,opt and 
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Fig. 22. Relation between visual extinction Ay and SFR 
for the s-MbN3Rs and the sb-MbN3Rs. Symbols are the 
same as in figure 18. Solid lines represent the fitted results 
with equation (22) for M* = IO^Mq at z = 1.5,1.0, and 0.6 
from top to bottom. Top-Left: For Ay-jR and SFRjrj^ijv 
of sb-MbN3Rs. Top- Right: For Ay -opt and SFRuv.corr of 
sb-MbN3Rs. Bottom-Left: The same as the Top-Left for 
s-MbN3Rs. Bottom-Right: The same as the Top-Right for 
s-MbN3Rs. 



8.1. Correlation among SFR, extinction, redshift, and 
mass 

Figures 22, 23, and 24 show the correlation between 
Ay and SFR for s-MbN3Rs, sb-MmN3Rs, MfNSRs, and 
MfBBGs. For all the subsamples, we have fitted correla- 
tions among Ay, SFR, z and M* with a form: 



SFR 
SFRo 



10"^''(l + z)'= 



10"M< 







(22) 



where SFRq is a normalization parameter for the SFR in 
the unit of Moyr"^. Their best-fitting planes arc deter- 
mined in the space of Ay, log(l + z), logM*, and \ogSFR 
as summarized in table 5 and shown in figures 22, 23, 



Fig. 24. Relation between visual extinction Ay and 
SFRuv.corr for the MfN3Rs and the MfBBGs. Solid 
lines represent the fitted results with equation (22) for 
M* = 10^" Mq at 2 = 1.5,1.0, and 0.6 from top to bottom. 
Left: For the MfNSRs. Right: For the MfBBGs. 

and 24. For the s-MbN3Rs, we can estimate both classi- 
cal visual extinction Ay with the SED fitting and calori- 
metric visual extinction Ay.jn converted from the calori- 
metric UV extinction Ajjy-iR with equation (19) for the 
best-fitted reddening curve in the optical-NIR SED fitting. 

We should remember that SFRm+jjv and Ay-jn are 
estimated from their MIR SED fittings while SFRjjy-corr 
and Ay-opt arc from their optical SED fittings. This de- 
notes that the trends of SFRjji^uy and Ay-jn arc inde- 
pendent of those of SFRuv;corr and Ay-^opt- As shown in 
figure 22, however, the correlation between SFRm+uv 
and Ay.j]:j is similar to that between SFRjjy.corr and 
Ay-^opt- Furthermore, the trends, among S F Rjjy .Qorrj 
■^V;opt and Zpiiotzy except M^, are similar to each other 
not only in the MbN3Rs and the MmN3Rs but also in 
all the subsamples including the MfN3Rs and MfBBGs as 
shown in figure 22 (see also table 5). All of them show 
the same trends; 1) larger SFR, larger extinction Ay, and 
2) lower redshift, larger extinction Ay. The second trend, 
increasing Ay with the cosmic time, suggests that increas- 
ing dust has been induced with chemical enrichment in 
evolving galaxies. 

8.2. Mass dependence and evolution of SFR 

Figures 25 and 26 show SFRir+jjv to Af* for the 
sb- and s/a-MbN3Rs and sb-MmN3Rs. We can see 
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Table 5. Fitting results with equation (22). 
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/3 
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sb-MbN3R 


0.4 


0.16 


4.2 


0.17 




±0.1 


±0.02 


±0.2 


±0.04 


s-MbN3R 


0.4 


0.19 


3.7 


0.13 




±0.1 


±0.02 


±0.2 


±0.04 


sb-MmN3R 


0.4 


0.05 


4.2 


0.16 




±0.1 


±0.02 


±0.2 


±0.03 


Ay -opt 










sb-MbN3R 


-1.1 


0.67 


5.2 


-0.19 




±0.2 


±0.05 


±0.5 


±0.08 


s-MbN3R 


-1.0 


0.72 


4.5 


-0.21 




±0.1 


±0.04 


±0.4 


±0.06 


sb-MinN3R 


-1.9 


0.9 


5.5 


-0.42 




±0.3 


±0.1 


±0.8 


±0.11 


MfN3Rs 


-3.3 


1.3 


6.5 


-0.71 




±0.2 


±0.1 


±0.4 


±0.06 


MfflBGs 


-2.3 


1.1 


5.6 


-0.19 




±0.2 


±0.1 


±0.1 


±0.06 



Values are means ± standard deviations. 
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Fig. 25. Stellar mass M* vs. SFRm^uy derived with the 
TIR and observed UV luminosities. The lines are the same as 
in figure 27 (see the text in detail). Symbols are the same as in 
figure 18. Left: For the sb-MbN3Rs. Right: The same as on 
the left for the star- forming components in the s/a-MbN3Rs. 



correlations between SFRmj^uv and M*, which are 
more obvious for the sb-MbN3Rs than for the s/a- 
MbNSRs. The correlation between SFRm+uv and Af* 
in the sb-MbN3Rs has also an evolutionary trend. The 
SFRinj^uy-Mjf correlation in the s/a-MbN3Rs is more 
scattered to the lower SFR side from that in the sb- 
MbN3Rs, in which the decrement becomes more obvi- 
ous in lower redshifts z < 0.8 even though the highest 
mass group around lO^^M© at higher redshifts z > 0.8 
tends to mostly follow that in the sb-MbN3Rs. They in- 
dicate that star formation activities in the s/a-MbN3Rs 
were quenched at z < 0.8 compared with those in the sb- 
MbNSRs. The indicated connection between quenching 
star formation and AGN activity around z ~ 0.8 will be 
discussed with the evolution of the stellar populations in 
section 9. 

The evolutionary trend and mass dependence on the 
SFR can be seen more clearly in figure 27, for all the 
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Fig. 26. The same as figure 25 for the sb-MmN3Rs. Left: 
For the sb-pP-MmN3Rs. Right: For the sb-fP-MmN3Rs. 



star-forming galaxies including the s-MbN3Rs, the sb- 
MmN3Rs, the star-forming MfN3Rs and the MfBBGs, 
which are plotted from the top to the bottom in red- 
shift intervals of 0.4 < z < 0.8,0.8 < z < 1.2, and z > 
1.2. The star-forming MfN3Rs and MfBBGs are selected 
as the population in the region of "blue cloud" in the 
([/ — V^) — A[/y color as remarked in subsection 6.4, which 
excludes extremely young objects fitted with the BC03 
models of age f < 0.1 Gyr. 

The mean SFR. for all the star-forming galaxies in- 
cluding both the s-MbN3Rs, sb-MmN3Rs, MfN3Rs, and 
MfBBGs is presented as an approximate function of A/* 
and z: 



SFR{MQyr-^\ :^ SFRu 
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(23) 



where the best-fitting plane in the space of logM*, log(l + 
z), and \ogSFR is determined as SFRh/Mqut-^ = 7.8, 
a = 0.74 ±0.01, and /3 = 2.7 ±0.1. In figures 25, 26, and 
27, the solid lines at z = 0.65,1.0, and 1.5 represent the 
best-fitting result. Roughly taking a ~ 1 and /3 ~ 3 for 
simplification, the growth time scale of the star-forming 
galaxies is typically derived as 

J^[yr]^10^"(l + z)-^ (24) 

which means the inverse of the sSFR. These SFR evo- 
lutionary properties and mass dependence are consistent 
with previous work for star-forming galaxies selected with 
general schemes (Noeske et al. 2007; Santini et al. 2009; 
Dunne et al. 2009; Pannella et al. 2009). 

9. Evolution of stellar populations 

9.1. Rest-frame optical color distributions 

The CMD is a powerful tool for studying galaxy evo- 
lution since it shows the bimodal galaxy distribution as 
early types that concentrate in a tight red sequence while 
late types distribute in a blue dispersed cloud, which have 
been studied not only in the local Universe with the SDSS 
(Blanton et al. 2003) but also in the distant Universe with 
the surveys for z < 1 (Bell et al. 2004; Faber et al. 2007; 
Borch ct al. 2006) and for z > 1 (Pannella et al. 2009; 
Brammer et al. 2009). Spitzer results have also shown 
that rest-frame colors of MlPS-detectcd objects, which 
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Fig. 27. Stellar mass M« vs. SFR for all the star-form- 
ing galaxies; the sb-MmN3Rs, MfN3Rs, and MfBBGs. The 
lines represent the fitting result (see the text in detail). 
Symbols are the same as in figure 6, except for the large dark 
and small thin ones represent the MIR-detected star-form- 
ing N3Rs (s-MbN3Rs and sb-MmN3Rs) and the MIR faint 
star-forming galaxies (star- forming MfNSRs and MfBBGs), 
respectively. Top: For z ^ 1.2. Middle: For 0.8 < 2 < 1.2 
Bottom: For 0.4 < 2 < 0.8. 



may overlap with the MbNSRs and the MmNSRs, fall be- 
tween the red sequence and the blue cloud up to z = 2 
(Cowie & Barger 2008; Salim et al. 2009; Brammer et al. 
2009). 

Thus, we also study rest- frame U — V colors for MbNSRs 
and MmNSRs in the field from their magnitudes interpo- 
lated with equation (4) . Figure 28 shows the U — V color 
without extinction correction for the stellar mass M^, for 
the sb-, s/a-, and agn- MbNSRs from top to bottom, where 
the left, middle, and right figures correspond to the red- 
shift intervals of 0.4 < z < 0.8, 0.8 < z < 1.2, and z > 1.2, 
respectively, and overplotted contours represent the rest- 
frame V detected galaxies, which are also shown with the 
MfNSRs in figure 29. 

The well-known "red and dead" population in a red se- 



Fig. 28. Stellar mass and rcst-framo color diagram of Af* vs. 
U-V tor the MbN3Rs excluded BL AGNs and EBOs. The 
color is not corrected for dust extinction. Typical errors of 
the mass and color are represented with crosses at the bottom. 
Large open and filled symbols and filled symbols in circles rep- 
resent s-,p-, and qp-BBGs, respectively, which arc the same 
as in figure 6. The plots are for 0.4 < 2 < 0.8, 0.8 < 2 < 1.2, and 
2 ^ 1.2 from left to right. Red and blue contours represent 
the rest-frame V detected galaxies with sSFR < O.lGyr"^ 
and sSFR > O.lGyr"^, respectively. The dotted line corre- 
sponds to Bell's critical line dividing the photometric samples 
into the blue cloud and the red sequence. Top: Largo colored 
and small gray symbols represent the sb-MbN3Rs and the 
2'-dctectcd galaxies, respectively. Middle: The same as the 
Top for the s/a-MbN3Rs. Bottom: The same as the Top for 
the agn-MbN3Rs. 

quence has a typical color of C/ — F ~ 2, which can be the 
reddest ones even in the local Universe with low Ay val- 
ues. However, redder populations ofU — V >2 can clearly 
be seen in all the diagrams not only for the MbNSRs but 
also for the MmNSRs and even for the MfNSRs as shown 
in figures 28, SO and 29. It suggests that not only the 
MbNSRs but also the MmNSRs and most of the MfNSRs 
suffer from dust extinction, which is already expected from 
the results in section 8. Thus, we have taken into account 
color corrections for these diagrams. Figures SI, S2, and 
SS represent intrinsic stellar color {U~V) — Aj/y corrected 
for extinction using Auv;opt derived from the optical-NIR 
SED fitting. In the top diagram of figure SI, we can see 
that the bimodal color distribution of the MbNSRs per- 
sists even up to z ~ 1.2, which is also consistent with the 
BBG classification as p-BBGs mainly appear in the region 
of a red sequence. As discussed in appendix l.S, unfortu- 
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Fig. 29. The same as figure 28 for the the rest-frame V 
detected galaxies and the MfNSRs. The red and blue small 
dots with contour represent the rest-frame V detected galaxies 
with sSFR< O.lGyr-l and sS FR > O.lGyr-^ , respectively. 
Top: For the rest-frame V detected galaxies. Bottom: For 
the MfNSRs. 




Fig. 30. The same as figure 28 for the MmNSRs. Top: Large 
colored symbols represent the sb-MmN3Rs. Bottom: The 
same as the Top for the nsb-MmN3Rs. 

nately, our Kg photometry is not deep enough to obtain 
a complete sample up to z> 1.2, which may be one of the 
reasons why the bimodal feature at z > 1.2 is not so clear 
compared with the results of Brammer et al. (2009). 

Figure 35 also shows the intrinsic colors of stellar com- 
ponents in the sb- and s/a-MbN3Rs, in which the col- 
ors are corrected with AuvjR- Few of them remain 
in the red sequence region, which comes from the ef- 
fect of the extinctions as discussed in subsection 7.3: 
{Auv-jr) ^ (Auviopt) + (AAuv) with (AAuv) 0.5 as 
long as fsK = 1- Furthermore, we cannot clearly see the 
bimodality in figure 35 compared with figure 31. The same 
trends can be also confirmed for the MmN3Rs as shown 
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Fig. 31. The same as figure 28 for color U — V — ^uv cor- 
rected using the extinction Aj/y.^pt, which is derived from 
the optical-NIR SED fitting. The' sb-,s/a-, and agn-MbN3Rs 
are plotted from the top to the bottom. 
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Fig. 32. The same as figure 29 for color U — V — Aj/i/ cor- 
rected using the extinction AijV;opt- 

in figures 31 and 36. Even if we have {/S.Auv) ^ by 
taking the proportional factor as fsk — 0.3, we just obtain 
a bluer mean color by adding ~ 0.5 without recovering 
the bimodality. It suggests that it is appropriate to cor- 
rect the color of a whole galaxy in the classical way with 
Auv;opt even for dusty systems as the MbN3Rs and the 
MmN3Rs. Thus, we will discuss the evolution of stellar 
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Fig. 33. The same as figure 30 for color IJ — V — Ayy cor- 
rected using tlie extinction A(7V;opt- 



Fig. 35. Tile same as figure 31 for color U — V — I^vv cor- 
rected using the extinction Auy-^jji^ which is derived with 
equation (20). Top: For the sb-MbN3Rs. Bottom: For the 
s/a-MbN3Rs, its remaining dusty starburst component can 
be extracted and then used to derive SFRm^jjy by fitting 
its MIR SED with the S&K models. 




Fig. 34. Histograms of the corrected color distribution with 
Ay derived from optical SED fittings. Height of histograms 
is normalized with the maximum number, which is displayed 
on the upper-left side. Top: For the rest-frame V detected 
galaxies, which are classified into two groups with sSFR 
larger and less than O.lGyr"^, which corresponds to a color 
U-V-A = 1.2. Bottom: For the sb-,s/a-, and agn-MbN3Rs. 

populations along with the results corrected with AiJV;opt 
as figure 31 in the following. 

We can see the evolutionary trends for sb-, s/a-, and 
agn-MbN3Rs in figures 31 and 34 and will summarize as 
follows. In all redshifts, the sb-MbN3Rs mainly stay in 
the area U — V — A <1. It can be interpreted as they are 
on-going in star formation. On the other hand, the agn- 
MbN3Rs start to appear even around U — V — A<la.tz> 
1.2 with large scatters in a mass range > 3 x lO^^M©, and 
relatively concentrate around U — V — A^^ 1 at 0.8 <z< 1.2 
as shown in the bottom diagrams of figure 31. The trend, 
as starbursts are bluer than AGNs at0.8<z<1.2, can be 
seen not only for the MbN3Rs but also for the MmN3Rs 
in figure 33. Figures 31 and 33 also show that the sb- 
MmN3Rs are relatively redder than the sb-MbN3Rs in the 
same stellar mass. It is reasonable as the former with less 
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Fig. 36. The same as figure 35 for the MmN3Rs. Top: For 
the sb-MmN3Rs. Bottom: For the non sb-MmN3Rs. 

luminous IR emissions may have smaller SFRs and quench 
their SF activities more effectively than the latter. 

As shown in the middle of figure 31, the s/a-MbN3Rs 
have a similar stellar population of the sb-MbN3Rs at 
z > 0.8 even though some of them were U — V ~ A > 1. 
Thus, a dusty star formation cannot suffer from the ob- 
scured AGN activity at z > 0.8 at least. On the other 
hand, at z < 0.8, it is similar to the agn-MbN3Rs at 
around [/ — y — A ^ 1. Passive populations such as the 
p-BBGs and on the red sequence are relatively major 
in s/a-MbN3Rs at 0.8 < z < 1.2. This suggests the co- 
existence of the obscured AGN with passive stellar popu- 
lations in their hosts, which are older than those in the sb- 
MbNSRs. Thus, the stellar populations become already 
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passive in massive galaxies such as the s/a-MbN3Rs of 
~ 1O^^M0 with obscured AGNs even at z ~ 1, and the 
obscured AGNs arc activated in less-massive galaxies of 
< IO^^Mq mainly with star formation at z < 0.8. These 
indications from the CMDs are consistent with the indica- 
tion about quenching SFRs of the s/a-MbN3Rs at z < 0.8 
from figures 25 and 26. It suggests that the evolution of 
AGN activity is more rapid than that of star-formation 
around z = 1 — 2 and that the co-cvolutionary link be- 
tween AGN activity and star formation in galaxies can be 
critical around z ~ 1. This seems to remind us of a sce- 
nario where AGN activity may quench the star formation 
in the distant universe to explain the downsizing. 

10. Evolution of AGN activities 

10.1. AGN activity traced at rest-frame 5 fim 
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Fig. 37. Correlation between stellar mass A^* vs. i/Lu{5 /im) 
at rest-frame 5 fim for the MbN3Rs. Symbols represent the 
same as in figure 6 except for the colors. The blue, green, and 
red symbols represent objects at the redshifts of 0.4 < z < 0.8, 
0.8 < z < 1.2, and z > 1.2, respectively. Diamonds repre- 
sent the spectroscopic BL AGNs. Dark and faint symbols 
represent objects with reduced < 0-5 ^^'^ > 0-5 ™ the 
MIR SED fitting, respectively. Top; For the agn-MbN3Rs. 
Bottom- Left: For the s/a-MbN3Rs. Bottom-Right: For the 
sb-MbN3Rs. 
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Fig. 38. The same as figure 37 for the MmN3Rs. Left: For 
the non sb-MmN3Rs. Right: For the sb-MmN3Rs. 

Figure 37 shows the rest-frame absolute monochromatic 
luminosities vL^-^^ at 5 /im to their stellar mass derived 
with equation (4), for the agn-MbN3Rs, the sb-MbN3Rs 
and the s/a-MbN3Rs. We can see an evolutionary trend 
as the agn-MbN3Rs at higher redshifts are more luminous 
at 5 /im. The mass and evolutionary trend for vL^-c, is 
roughly summarized as: 
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(25) 



In figure 37, the data points of the MbN3Rs arc plotted 
in rcdshift intervals of 0.4 < z < 0.8,0.8 <z< 1.2, 1.2 < z, 
respectively. The solid line for agn- and s/a-MbN3Rs at 
z = 0.6, 1.0, and 1.5 represents the best-fitted plane of 
log(Lii/L0) = 8.9±O.l, a = 0.15 ±0.08, and /? = 6.3 ± 0.4 
in equation (25). On the other hand, the solid lines for sb- 
MbN3Rs represent the plane of a = 0.74 and (3 = 2.7, which 
are the same expornents of equation (23) for star-forming 
activities in subsection 7.1. Thus, the 5-/im emissions 
of the sb-MbN3Rs are dominant from their dusty star- 
forming regions. 

The trends for the 5-/im luminosities do not depend 
on their selections since the i^L^-^ of the MbN3Rs is 
larger than the limiting luminosities for their detection 
as ~ 10^ 6 X 10^ and 2 x IO^^Lq at z = 0.5, 1.0, and 1.5, 
respectively. Thus, the emission process at the rest-frame 
5 fim for the agn- and the s/a-MbN3Rs should be differ- 
ent from the dusty star-forming origin for the sb-MbN3Rs. 
This supports the supposition that MIR emissions from 
the agn-MbN3Rs are dominated by obscured AGN. Thus, 
it is expected that the luminosity vL^y-^ may dom- 

inate in the bolometric luminosity of the obscured AGN 
Lboi for the s/a- and agn-MbN3Rs. 

The evolutionary trend of I'L^.^ in the agn-MbN3Rs is 
more rapid than that of SFRs in the sb-MbN3Rs. The 
vL^;5 in the agn-MbN3Rs at 0.8 < z < 1.2 and 0.4 < 
z < 0.8 also becomes less scattered than those at z > 1.2, 
which can be attributable to quenching the AGN activities 
around z ~ 1 . 

The analysis with vL^.^^ can also apply for the MmN3Rs 
as shown in figure 38 even though they are biased since 
most of them are detected in critical threshold at the rest- 
frame 5 /im. The sb-MmN3Rs and the non sb-MmN3Rs 
show roughly a similar trend to those of the sb-MbN3Rs 
and the s/a-MbN3Rs compared with the agn-MbN3Rs, 
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which suggests that AGN activities are not dominant in 
the MmNSRs. 

10.2. Super Massive Black Holes in AGNs 

Consider accretion onto the Super Massive Black Hole 
(SMBH) with a mass Mbh for an AGN. The bolomctric 
luminosity L^oi of the emission from the AGN is charac- 
terized by an Eddington ratio A as L^oi = AL Edd and the 
Eddington luminosity Lsdd- 



LEdd [erg s'^] = 1.26 x 10 
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Introducing a luminosity ratio = i/L^.r^/ L{,oi, we can 
obtain a constraint for the mass of the SMBH from the 
5-/im luminosity i^L^-^ as 

Mbh 



3.06 X IO6M0 



■'■5 iniir 



(27) 



We can see vL^-^^ ~ lO^^Lo, and 3 x \Q^^Lq at maximum, 
and minimum, respectively, for the agn-MbN3Rs of ~ 
lO"A/0 at z > 1.2 in figure 37. It means that SMBHs 
had already grown in a mass range of ~ 1O®A~^/5~^M0 — 
3 x WX-^J^^Mq > IO^Mq in a host of A/, ~ 10"Mq 
selected as the agn-MbN3Rs at z > 1.2. Furthermore, 
as long as A/5 > 10""'^, the BH mass in a maximum 
luminous agn-MbN3R of ^ IO^Az/q at z > 1.2 does not 
exceed 3 x 10^ Af©, which is a similar value attributable to 
the relation between the mass of SMBH and the stellar 
mass in the local spheroid as derived by Haring & Rix 
(2004): 



log Afsff = -4 . 1 2 + 1 . 1 2 (log 



(28) 



If the mass relation between SMBH and host was already 
established at z > 1.2 as this equation (28) in the local 
Universe, then, the AGN luminosity in the agn-MbN3Rs 
at z > 1.2 could be near the Eddington limit A ~ 1, which is 
consistent with consideration about SMBH growth history 
as follows. 

Introducing an energy conversion efficiency e = 
Lboil{M(?) for mass accretion with M onto an SMBH, 
the 5-/im luminosity vLi,.^ also implies a constraint for a 
mass growth ratio Mbh of SMBHs as 
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BH 
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(29) 



where we have used Mbh = (1 ^ ^)M and vL^-^^ = f^Ltoi. 
The thin-disk accretion model has predicted the efficiency 
e should be in the range ~ 0.04 — 0.31, depending on the 
spin of SMBH. The growth time scale of SMBHs is also 
represented as 

:^^[yr]=4.5xl0*^A-i^^. (30) 
Mbh^ 

Then, the SMBH growth time could be less than the time 
of ~ 2.5 Gyr from z=2 to z=l for A/e > 0.2, which can 
be satisfied as long as it is near the Eddington limit with 
A > 0.06 for assuming the thin-disk accretion model with 
e — 0.04 — 0.31. Thus, it is consistent with the suggestion 



that SMBHs had already grown in the maximum 5-//m 
luminous agn-MbN3Rs of Af* ^ 10" A/© at z > 1.2 and 
their AGN luminosity could be near the Eddington limit 
A ~ 1 as discussed above. 

10.3. AGN activity and star formation 

Even though the sb-MbN3Rs and the agn-MbN3Rs are 
different populations, we could obtain simple formulae 
representing the activity of star formation and obscured 
AGN in the MbN3Rs up to z ~ 2 with functions of their 
stellar mass and redshift as equations (23) and (25). The 
5-/xm luminosity of the agn-MbN3Rs is rapidly quenching 
as equation (25) compared with that of sb-MbN3Rs as 
equation (23). It means that the evolutionary time scale 
of AGN activity is shorter than that of star formation for 
the MbN3Rs selected with the AKARI MIR photometry. 

The rapid evolutionary trend of AGNs in the MbN3Rs 
was already expected from its color of stellar population as 
discussed in section 9. Thus, this evolutionary trend in the 
time scales, which distinguishes between the agn-MbN3Rs 
and the sb-MbN3Rs, is consistent with the evolutionary 
feature of the MbN3Rs around z ~ 1 in their distribution 
on the color- mass diagram as shown in figure 35. 

11. Discussions 





8 






ff 
















a 


P 
















10 20 20 50 1 00 200 500 
wavelength (/im) 



10^ 10'" 



Fig. 39. Left: Averaged rest frame IR-Submm SEDs for 
the sb-MbN3Rs. Red filled circles, green open squares, 
and blue open circles represent monochromatic luminos- 
ity uLi, averaged for the sb-MbN3Rs which fit well with 
large, medium, and small optical depth S&K models of 
^V SK = 17.9,6.7, and 2.2, respectively. Red solid, green 
dotted, and blue dashed lines represent model SEDs with 
^V;SK = 17.9,6.7, and 2.2, respectively. The model SEDs 
have been convoluted within a 3 fim wavelength corre- 
sponding to a mean width of the AKARI/IRC bands. 
Right: Expected monochromatic luminosity at rest-frame 
600 fim; !^L^(600 ^m) vs. monochromatic luminosity ratio 
at rest-frame 10 fim to 7.7 /im; uLi,{lO fim)/uL,y{7.7 fim) for 
a case with TIR luminosity Lm = 5 X lO'^^ Lq. Black filled 
circles, open squares, and open circles represent monochro- 
matic luminosities i'L^(10 fira) /vLviJ .7 fim) from the con- 
voluted S&K models with Ay-sK = 17.9,6.7, and 2.2, re- 
spectively. We have taken other S&K model parameters as 
L<°t = 10i° l,10li l,andl0l2 l Lq from right to left and fixed 
ios/i*°* = 0.6. 



11.1. FIR and Submillimeter Follow-ups 

As discussed in subsection 7.3, the s-MbN3Rs and 
the sb-MmN3Rs with a lower ratio of monochro- 
matic luminosity at rest-frame 10 /zm to that at 



26 



Hanami et al. 



[Vol. , 





MbN3H 

0.4 S z < O.B 




MinN3R Qsb •' 

4 < z < O.B Anbn sb-' 




Fig. 40. BzK color diagram for tlie MbNSRs and MmNSRs. 
Large darlc and faint symbols represent galaxies at 2 > 1.2 
and z < 1.2), resepectively. The small dots represent the MIR 
faint galaxies with K < 21. Left: For the MbN3Rs. Red 
open circles, green open squares, blue open triangles, and 
magenta filled diamonds represent sb-, s/a-, agn-MbN3R, and 
BL-AGN, respectively. Right: For the MmNSRs. Rod open 
circles, green open squares and magenta filled diamonds rep- 
resent sb-, non sb-MmN3R, and BL-AGN, respectively. 
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Fig. 42. The same as figure 40 for the uVi color diagram. 
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Fig. 41. The same as figure 40 for the uRJ color dia- 
gram. Large dark and faint symbols represent galaxies at 
0.8 < z < 1.2 and z < 0.8, respectively. Small dots represents 
the MIR faint galaxies writh J < 22. 

7.7 fiia; z/Ljy(10 ^m) / vL^{7 .7 /im) tend to have 
larger extinctions assuming denser dust distributions 
(see in figure 21). Larger extinctions, indicated with 
vLy{lQ ixm)/vL,j{7.7 /im), can make reprocessed photon 
emissions dominant at longer wavelengths. This trend can 
be confirmed in the left plot of figure 39, which shows 
the average IR/Submm SEDs for the s-MbN3Rs fitted 
well with large, medium, and small optical depth param- 
eters in S&K model; Ay-SK = 17.9, 6.7, and 2.2, respec- 
tively. The s-MbN3Rs with a larger optical depth tend 
to show not only a deeper dip at 10 ^va but also a larger 
flux in FIR and submillimeter wavelengths. This feature 
can be confirmed directly with observations with FIR de- 
vices such as the Photodctector Array Camera (PACS) 
and Photometric Imaging Receiver (SPIRE) on board 
Herschel and submillimeter telescopes as MAMBO/IRAM 
and ALMA. Figure 39 shows the expected monochro- 
matic luminosity at rest- frame 600 ^m; i^Li,(6GG /im) vs. 
i/Ljy(10 ^Lva) / vL,j[l .7 /im) as the trend between extinc- 
tions and FIR/Submm wavelength excess. 

Observations with ground-based submillimeter tele- 
scopes as SCUBA and MAMBO have reported that 
so-called high-z SMGs tend to have SEDs reproduced 
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Fig. 43. Color distribution of the MbNSRs along the strip 
between the dashed lines in figures 40 and 41. Top solid line 
histogram indicates the number of galaxies detected in the 
NIR band. Red slashed, filled, and outline histograms repre- 
sent agn-, s/a-, and sb-MbN3Rs, respectively. Blue slashed 
and outline histograms represent and non sb-MmNSRs and 
sb-MmNSRs, respectively. 

with lower dust temperatures when compared with local 
ULIRGs (e.g. Chapman et al. (2004)), which suggests 
that detection of the SMGs with the ground-based submil- 
limeter telescopes could be biased and there are missing 
populations of star-forming galaxies with hotter dust as 
Optically Faint Radio Galaxies (OFRGs). Recently the 
PACS and SPIRE/77ersc/ie/ observations have revealed 
that ULIRGs at z ~ 2 detected with the MIPS / Spitzer 
consist of both cold and warm populations correspond- 
ing to the SMGs and the OFRGs (Magdis et al. (2010); 
Elbaz et al. 2011), respectively. Even though our sam- 
ples of MbN3Rs and MmN3Rs may not overlap with the 
Herschel detected ULIRG samples since our samples were 
detected in the MIR bands and correspond to LIRGs at 
z ~ 1, the variety of vL^{\0 ^m)/h'L^{7.7 fim) suggests 
that LIRGs at z ~ 1 may consist of various populations 
with cold and warm dusts similar to the ULIRGs at z ~ 2 
such as the SMGs. The forthcoming and future obser- 
vations with Herschel and submillimeter telescopes will 
confirm these features as predicted in figure 39. 

Adding these photometric data at these longer wave- 
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lengths for the MbNSRs can also serve as a strong con- 
straint in the SED analysis with direct estimation of their 
TIR luminosities, which can overcome an ambiguity left in 
deriving absolute TIR luminosity by introducing the pro- 
portional factor fsx in section 7. Furthermore, this makes 
a distinction between starburst dominant and AGN domi- 
nant with the FIR/Submm luminosity with the MIR one, 
which is an alternative way to determine the MIR color 
criteria for selecting s- and agn-MbN3Rs in this work. 
The alternative selection for AGNs can be also effective 
to study the census of AGN populations. As shown in 
subsection 6.1, the CXO found that some agn- and s/a- 
MbNSRs are associated with X-ray sources. On the other 
hand, there are still some populations in the agn- and 
s/a-MbN3Rs escaping from the X-ray detection with the 
CXO. This suggests that obscured AGNs in the agn- and 
s/a-MbN3Rs are classified to two subpopulations associ- 
ated with and without X-ray sources. The former and the 
latter are possible candidates of obscured AGNs in the 
Compton-thin range of TV// < 10^*cm^^ and the Compton- 
thick (CT) range of Nh > lO^'^cm"^, respectively. Thus, 
the comparison of MIR-FIR/Submm SEDs between the 
CT and non-CT AGN candidates is expected to supply 
constraints for physical states around the SMBHs in the 
obscured AGNs. 

11.2. Optical colors of MbNSRs and MmNSRs 

Most of the MbNSRs and MmNSRs are also classi- 
fied as BBGs. Then, we have also compared the sb- 
, s/a-, and agn- MbNSRs with their color on the two- 
color diagrams related to the BBG criteria derived in 
appendix 1. Figures 40, 41, and 42 show that the s- 
MbNSRs and the sb-MmNSRs mainly concentrate around 
{B~z,z~K) :^ (2,2.5), (u- R,R- J) ~ (2,2.5), and 
{u-V,V-i):^ (1.5, 1.5), while the agn-MbNSRs are rela- 
tively more dispersed than the s-MbN3Rs. These BfiTZ 
colors of the s-MbNSRs and the sb-MmN3Rs are con- 
sistent with those expected from the reddenings with 
heavy extinctions as shown by Daddi et al. (2004) and 
in appendix 1. We detected ~ 60 s-BzKs and ~ 100 
s-uRJs with K < 21.5 and J < 22 in the color region 
within A{B - K) ~ ±0.5 and A(w - J) ±0.5 around 
[B - z,z- K) ~ (2, 2.5) and [u- R,R- .J) ~ (2, 2.5), re- 
spectively, where ~ 15 and ^ 40 out of these s-BzKs (s- 
uRJs) are classified as s-MbNSRs as shown on the left 
in figure 4S. Thus, roughly, we can detect dusty star- 
bursts at z 1 — 2 with a frequency one fifth of s-BzKs and 
one third of s-uRJ. For NIR bright BBGs with K < 20.5 
and J < 21, the fractions of MbNSRs and MmNSRs in 
BBGs increase as shown on the right in figure 4S. In 
particular, roughly half of the uRJs with J < 21 around 
(u — R, R — J) ~ (2, 2.5) are LIRGs around z ^ 1, which 
are selected in the NEP deep field. We should note that 
this NEP deep field is a blank field. Thus, these kinds 
of BBGs constrained with NIR brightness and colors in 
any blank field can be mostly detected as LIRGs around 
z ^ 1 — 2, which are interesting targets for a future un- 
biased survey for galaxies with on-going, up-coming, or 
future FIR/Submm telescopes as Herschel and ALMA. 



These future studies will be also essential to confirm the 
trends for the populations marginally detected with the 
AKARI as the MmNSRs. 

12. Conclusions 

We have studied BBGs/IRBGs in three redshift ranges 
of z = 0.4 -0.8, 0.8 -1.2, and > 1.2 with MIR muhiband 
photometry, which revealed their star- forming/ AGN ac- 
tivities with/ without PAH emissions deriving SFR from 
their TIR luminosities. 

• The MIR colors and SED analysis can classify 
MIR bright NS Red galaxies (MbNSRs) into three 
classes; starburst dominates (sb-MbNSRs), star- 
burst /AGN co-existence(s/a-MbNSRs), and AGN 
dominates (agn-MbNSRs). The ratio of rest- frame 
luminosities at 5 and 10 /im to that at 7.7 /im is sys- 
tematically higher in the s/a-(agn-)MbNSRs than in 
the sb-(s/a-)MbNSRs, which is consistent with the 
picture of the s/a-MbNSRs as the starburst / AGN 
co-existent subpopulation in the dusty star-forming 
MbNSRs (s-MbNSRs). 

• MIR marginally-detected NSRs (MmNSRs) can be 
also classified into two classes; sb-MmNSRs as mim- 
ics of sb- MbNSRs, and non sb-MmNSRs. 

• Their rest-frame 7.7-/xm luminosity is a good tracer 
of TIR luminosity as the PAH emission dominates 
for the sb- and s/a-MbNSRs and sb-MmNSRs which 
are classified as dusty star-forming galaxies even up 
to z = 2. 

• The SFRjji+ijv derived from the total bolomet- 
ric luminosity Ljfi^jjv of the sb-MbNSRs/MmNSRs 
shows a correlation that is nearly proportional to 
the stellar mass M*, which is consistent with that 
of the extinction corrected UV star formation ratio 

SFRiJV;corr- 

• The SFRiR+uv of the s/a- MbNSRs tends to be 
smaller than that of the sb- MbNSRs at z < 0.8. 

• The IR-derived specific SFR of the sb- 
MbNSRs/MmNSRs rises with the redshift at 
all stellar masses. 

• The ratio of the rest-frame luminosity at 10 fim to 
that at 7.7 fim of the sb-MbNSRs/MmNSRs may 
trace the optical depth due to the dust extinction. 

• The AGN activity could be traced with the rest- 
frame 5-/im luminosity for the agn- and the 
s/a-MbNSRs. The evolutionary trend (1 + z)^ 
of the AGN luminosity for the agn-MbNSRs is 
more rapid than the (1 + z)^ of SFR for the sb- 
MbNSRs/MmNSRs. 

• SMBHs in the galaxies have grown to ~ S x 10^ Mq 
in maximum luminous populations with lO^^L© and 
10"Mo of the MIR selected AGNs at z > 1.2, which 
suggests the mass relation between the SMBH and 
its host is already established at z ~ 1 — 2 and the 
same as that at the present day. 

This study is based on data mostly collected with 
the AKARI IR astronomical satellite and at the Subaru 
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Appendix 1. Classification with Balmer/4000 A 
breaks 

In order to identify galaxies at various redshift intervals, 
we have generalized the two-color technique for detect- 
ing the Balmer/4000 A break, which is a major feature 
in the SED of galaxies. The Balmer and 4000 A break 
are often treated as a single feature, owing to their simi- 
lar wavelengths and the overlapping populations of galax- 
ies having these features. However, the breaks originate 
due to different physical processes, and behave differently 
as a function of age and population. Both breaks result 
from absorption in stellar atmospheres. The 'classical' 
Balmer break at 3648 A marks the termination of the 
hydrogen Balmer series, and is strongest in A-type stars. 
This means that it is the most prominent feature for inter- 
mediate ages of 0.3 — 1.0 Gyr in a Single Stellar Population 
(SSP). The Balmer break strength, Db, is defined as the 
ratio of the flux density F,, in the 3500-3650 A and 3800- 
3950 A bands around the break (Balogh et al. 1999). On 
the other hand, the 4000 A break arises because of an ac- 
cumulation of absorption lines of mainly, ionized metals. 
Its strength is defined using an index, Z)„(4000), that is 
based on the continuum regions, and the ratio of the flux 
density in the 3850-3950 A and 4000-4100 A bands. As 
the opacity increases with decreasing stellar temperature, 
the 4000 A break strength, L'„(4000), increases for older 
ages, and is largest in old and metal-rich stellar popu- 
lations. However, the Balmer break strength Db does 
not monotonically increase with age, instead reaching a 
maximum at intermediate ages (Kriek et al. 2006). The 
metallicity has a more minor influence for ages less than 
1 Gyr (Bruzual & Chariot 2003). Thus, both break fea- 
tures around 4000 A include information about the stellar 
populations in a galaxy. Although these breaks are in- 
dividually unresolved from broad-band photometry, the 
degenerate Balmer/4000 A break still provides a robust 
feature in the SEDs with galaxies of an age t > 500 Myr. 
In the following subsections, we will demonstrate that the 
Balmer/4000 A break can be used to select galaxies at 
various redshift intervals of interest to this work. 
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Table 6. The Parameters for BBG Selection 



Bands 








^min 


•^max 


uVi 


-0.2 


1.5 


1.0 


0.4 


1.2 


uRJ 


0.0 


2.5 


1.0 


0.8 


1.8 


BzK 


-0.2 


2.5 


1.0 


1.4 


2.5 


BRN3 


0.7 


3.0 


3.0 


0.8 


1.8 


uRm 


0.0 


3.0 


2.0 


0.8 


1.8 


Bzm 


-0.2 


2.5 


1.4 


1.4 


2.5 


uzNi 


-0.5 


2.5 


1.2 


1.4 


2.5 




2 4 6 2 4 

B - z' B - z' 



Fig. 44. BzK diagram in the field. Symbols of cross bar at 
(B — z' , z' — Ks) = (3,4.5) represent typical errors of the color. 
Magenta diamonds represent BL AGNs, in which squared 
diamonds represent lower-redshifted ones at 2 < 1.4. Left: 
For the spectroscopic samples, red circles represent sources 
at Zspec = 1.4 — 2.5. Right: The same as on the Left for the 
photometric samples, where blue cross and yellow circles rep- 
resent the stellar objects selected with the scheme described 
in subsection 3.4, respectively, and the reference stars of the 
BPGS star atlas. Red circles and squares represent sources 
at Zphat = 1-4 — 2.5, and at Zpf^^t > 2.5 respectively. 

A. 1.1. The BzK, uRJ , and uVi color-color diagram 

Daddi et al. (2004) proposed that z > 1.4 star- forming 
BBGs (s-BBGs) should occupy the upper-left region in a 
two-color B — z and z — K (BzK) diagram, as the bluer 
and redder colors trace UV excess and the Balmer/4000 
A break feature, representing the star forming activity 
and the old stellar component, respectively. We could 
also show that the BzK technique is efficient for selecting 
galaxies at 1.4 < Zphot < 2.5, as illustrated in figure 44. 

We will argue that this kind of two-color technique can 
be extended to selecting s-BBGs at certain redshift inter- 
vals with other sets of three photometric observed broad 
bands, in which their wavelengths are nearly in an equal 
ratio as a geometric series. By introducing B, /i, and TZ as 
their AB magnitudes of general bluer, medium, and red- 
der bands, resepectively, and defining a two-color index 
B^TZ generalized in BzK: 

BfiTZ={n-TZ)-hBM'3~^^), (Al) 

it follows that s-BBGs at certain redshift intervals can be 
selected by a criterion: 

Biin>ht.n. (A2) 

and passively evolving BBGs (p-BBGs) can be selected 
with a criterion: 
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Fig. 45. The same as figure 44 for the uRJs. Symbols of 
cross bar at («* — R,R — J) = (3,4.5) represent typical errors 
of the color. Red circles and squares represent sources at 
^phot = 0.8 — 1.8 and at Zp^cjt > 1-8, respectively. Left: For the 
spectroscopic samples. Right: For the photometric samples. 
Only one fourth of all photometric sources detected in the 
uRJ bands are plotted to reduce crowding in the diagram. 




♦ BL AGN(z >0.4) 
ffiBL AGN(z;;^°<0,4) 




Fig. 46. The same as figure 44 for the uVis. Symbols of 
cross bar at (n* — V,V — i') = (1, 3) represent typical er- 
rors of the color. Red circle and square represent sources 
at Zspec = 0.4 — 1.2 and Zspec > 1-2, respectively. Left: For 
the spectroscopic sample. Right: For the photometric sam- 
ple. Only one thirtieth of all photometric sources detected in 
the uVi bands are plotted to reduce crowding in the diagram. 



(A3) 



Even though wc can obtain only the upper hmit in their B 
magnitudes, we can select quasi passively evolving BBGs 
(qp-BBGs) with a criterion: 



(A4) 



We confirm the validity of this generalized BBG tech- 
nique as shown in figures 45 and 46. These generalized 
criteria are efficient ways to select photometric s/p/qp- 
BBGs at 0.8 < z < 1.8 and 0.4 < z < 1.2 in the field for 
the uRJ and uVi filter sets. These particular color crite- 
ria also allow efficient selections of the BBGs not only for 
1.4 < z < 2.5, but also 0.8 < z < 1.8 and 0.4 < z < 1.2. 

In order to study the validity and physical meaning 
of the two-color criteria that are generalized from the 
phenomenologically established BzK technique of Daddi 
et al. (2004), we used the BC03 models to reproduce 
the blue B — n and red colors ji — TZoi model galaxies, 
which are described and shown in figures 57 and 58 in 
appendix 5.2. The result from the stellar population syn- 
thesis study is essentially similar to the cases of BzKs, as 
summarized below. 




Fig. 47. Photometric rcdshift distributions for the BBGs in 
the field. The white, shaded, and dark area represent the 
S-, qp-, and p-BBGs, respectively. Left: For the classical 
BBGs. Right: The same as on the left for the extended 
BBGs. Left-Top: Lower, upper solid line and dashed line his- 
togram represent the uVis, the z'-detected galaxies, and the 
rest-frame V-detected galaxies, respectively. Left-Middle: For 
theuRJs. Left-Bottom: For the BzKs. Right: For the uRN3s. 
Right-Upper-Middle: For the BRN3s. Right-Lower-Middle: 
For the uzNSs. Right-Bottom: For the BzNSs. 



In the redshift ranges for each BBG selection in ta- 
ble 6, the galaxies with ongoing star formation appear in 
the upper-left region B^TZ > /s^tj of the generalized BfiTZ 
color diagrams. The duration of the star-formation (age) 
has little infiuence on the blue color B — ^, while the redder 
color fi — TZ increases with age. Only very young starbursts 
with ages less than 10 Myr, and without underlying older 
stellar populations, are located around the lower-left re- 
gion, falling just below the line of B^TZ = fts^n- On the 
other hand, the colors of normal galaxies, obtained from 
the CWW SED templates, fall outside the color area of 
the BBGs in the redshift ranges for each of the BBG selec- 
tions. This means that the generalized two-color criteria 
can also be used to exclude the lower-z objects. 

The reddening direction lies approximately parallel to 
the line defined by Bj-iTZ — fsfiTZ, since the dust extinction 
can be approximately represented by a power-law function 
of wavelength, implying that the reddening is similar in 
the B — ^ and /i — 7?. colors for star-forming galaxies in the 
various redshift ranges. This means that the generalized 
criteria BfiTZ > fi3f_in for selecting BBGs are also robust 
from their dust reddening. 

Only a few galaxies are found as faint blue objects 
at K — ^ ~ 0, suggesting that purely un-reddened, star- 
forming galaxies are rare, not only for s-BzKs at z ~ 2, 
but also s-uRJs at z ~ 1. On the other hand, we could 
find some fraction of these nearly extinction-free popula- 
tions in s-uVis at z ~ 0.6. 

Figure 47 shows the distributions of photometric red- 
shifts Zphot for all the z'-detected galaxies and the BBGs. 
Total numbers of detected uRJs and uVis as the lower-z 
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Fig. 48. The same as figure 2 for classical BBGs. Typical er- 
rors of the magnitude and color are represented with crosses 
at the bottom. The vertical dotted and steep slope dotted 
lines represent the limiting absolute V and !7-band magni- 
tudes at the mean redshift, respectively. In general, colored 
open, solid, and open-solid symbols represent objects clas- 
sified as S-, p-, and qp-BBGs, respectively. Basically, red 
circles, green triangles, and blue stars represent objects se- 
lected with the two-color criteria as BzKs, uRJs, and uVis, 
respectively. Black-large open, solid, and open-solid symbols 
represent mean values for the s-, p-, and qp-BBGs, respec- 
tively. Crosses on the black symbols represent their stan- 
dard deviations. Only a part of samples arc plotted to re- 
duce crowding in the diagram, whose fractions are presented 
as numbers in the right-bottom region of all figures. Top: 
The BzKs at 0.4 < 2 < 1.2, 1.2 < 2 < 1.4, 1.4 < z < 2.5, and 
z > 2.5 are plotted from loft to right. For BzKs also clas- 
sified as uRJs at 1.2 < z < 1.4, their symbols are changed 
to green triangles from red circles. For BzKs also classified 
as uRJs at 0.4 < 2 < 1.2, the symbols are changed to blue 
stars from red circles. Middle: The uRJs at 0.4 < z < 0.8, 
0.8 < z < 1.4, 1.4 < z < 1.8, and 1.8 < z < 2.5 are plotted from 
left to right. For uRJs also classified as BzKs at 1.8 < z < 2.5, 
their symbols are changed to red circles from green triangles. 
For uRJs also classified as uVis at 0.4 < z < 0.8, their sym- 
bols are changed to blue stars from green triangles. Bottom: 
The uVis at 0.4 < z < 0.8, 0.8 < z < 1.2, 1.2 < z < 1.8, and 
1.8 < z < 2.5 are plotted from left to right. For uVis also 
classified as uRJs at 1.2 < z < 1.8, their symbols are changed 
to green triangles from blue stars. For uVis also classified as 
BzKs at 1.8 < z < 2.5, their symbols are changed to red circles 
from blue stars. 

BBGs are greater than that of detected BzKs while the 
fractions of low-z interlopers in the uRJs and uVis are 
in the same order as those in the BzKs as shown in fig- 
ure 47. Thus, all the BBG selections can be effective to 
exclude low-z interlopers at z < 0.4, z < 0.6, and z < 1.2 
with the uVi, uRJ, and BzK criteria. These selected-out 
redshift ranges are consistent with the expected values of 
z < 0.4, z < 0.8, and z < 1.2 from the stellar population 
models (see in appendix 5). While the uVis are selected 
out at z > 1.2 as expected with the model tracks (see 
in appendix 5), the uRJs and BzKs gradually decrease 
in population at z > 1.2, and > 1.6 lower than expected. 



Low efficiency for BBG selections at the high redshift side 
may be caused from the shallowness of our ground-based 
NIR photometry. The limiting magnitudes of J < 21 and 
Kg < 22 with the KPNO/FLMG are shallow compared 
with studies for the BzKs in other survey fields (Kong 
et al. 2006; Lane et al. 2007; Quadri et al. 2007; Blanc 
et al. 2008). Thus, we take into account this limitation of 
the ground-based NIR photometry in the dataset when we 
study galaxies at z ~ 0.8 — 2 as uRJs and BzKs. We try to 
overcome the limitation due to the shallow ground-based 
NIR photometry for the uRJs and BzKs. We further ex- 
tend the two-color criteria used for detecting the Balmer 
in combination with AKARI NIR photometry as discussed 
in appendix 1.3. 

A. 1.2. Rest frame color of the BBGs 

It is worth comparing the classification with the two- 
color criteria detecting the Balmer break for the BBGs 
and their rest frame colors. As shown in the top and 
middle of figure 48, most of the BzKs at 1.2 < Zphot < 1-4 
are classified as uRJs while only a few of the uRJs at 
1.4 < Zphot < 2.5 are classified as BzKs. This is related 
to the shallowness in our ground based NIR photometry, 
which is not deep enough to select the population around 
z ~ 2 as BzKs as remarked in the previous section. For 
our sample up to z ~ 2, the uRJ criteria are more useful 
than the BzK criteria. We can see the bimodality even 
for the BBGs again as the boundary between the blue 
cloud and the red sequence is consistent with that by Bell 
et al. (2004), which are plotted as nearly horizontal dotted 
lines in figure 48. Even though the bimodality of BBGs 
at z > 1.2 is not so clear compared with that at z < 1.2, 
the mean colors of the s-BBGs and the p- and qp-BBGs, 
represented as black open symbols and solid and open- 
solid symbols in figure 48, are in the blue cloud region 
and in the red sequence, respectively. Thus, the uRJ and 
uVi criteria can reasonably classify the the star forming 
and the passively evolving populations as the s-BBGs and 
the p- and qp-BBGs at z < 1.8 in the field at least. It 
confirms that the generalized BBG criteria are also useful 
in order to distinguish a star-forming population from a 
passively evolving one. 

A. 1.3. Salvaging BBGs from NSRs at z > 0.8 

As shown in appendix 1, the BBG techniques employing 
not only the original BzK color but also the uRJ color can 
select galaxies around the redshift desert. However, we 
could not completely select these BBGs from the ground- 
based photometric dataset in the field since the ground- 
based JKs photometry with the KPNO/FLMG is not 
deep enough to pick up the IRBGs detected with the 
AKARI as seen in table 2. Since BBGs and IRBGs are 
almost the same populations at z > 0.4 as discussed above, 
their statistics can confirm the limitation of the JK depth. 
Roughly one third of N4Rs and N34 bumpers could not 
be identified as BzKs or uRJs as seen in the comparison 
between the left plot in figure 47 and the top-left plot 
in figure 4. Even though the limiting magnitude in the 
N3 band is comparable to those in J and K photome- 
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Fig. 49. Two-color diagrams for the extended BBGs se- 
lected from the N3Rs, in which the band is taken as 
the reddest band TZ in the BBG criteria. Red crosses at 
{B — n, fi — TV) = (4,4.5) represent typical errors of the color. 
Only a part of the N3Rs are plotted to reduce crowding in 
the diagram, whose fractions are presented as numbers at 
{B — 11,11 — TV) = (5,4.5) in all figures. Diamonds represent BL 
AGNs. Large and small symbols represent the spectroscopic 
and the photometoric samples, respectively. Top-Left: The 
same as figure 44 for the BzNSs. Circles and squares represent 
objects aX z = 1.4 — 2.5 and z > 2.5, respectively. Top-Right: 
The same as figure 44 for the uzN3s. Circles and squares 
represent objects at z = 1.4 — 2.5 and z > 2.5, respectively. 
Bottom-Left: The same as figure 45 for the BRN3s. Circles 
and squares represent objects at 2 = 0.8 — 1.4 and z > 1.4, 
respectively. Bottom- Right: The same as figure 45 for the 
uRN3s. Circles and squares represent objects at z = 0.8 — 1.4 
and z > 1.4, respectively. 

try as seen in tabic 2, the N3 band is effective to de- 
tect emissions around the rest-frame 1.6-//m IR bump for 
the objects at 2 > 1 more than the J and K bands. In 
order to salvage these missing BBGs as non-BzK(non- 
uRJ) in N3Rs around the redshift desert, we used A^3 
photometry instead of the ground-based NIR ones for 
BBG selections with BzN3/uzN3 {uRN3/BRN3) col- 
ors as shown in figure 49. The numbers of salvaged 
BzN3s and uzN3s (uRN3s and BRN3s) were greater than 
those of BzKs(uRJs) even in this dataset as shown in fig- 
ure 47. Hereafter, we often treat BzN3s and uzN3s (uRN3 
and BRN3) as an extended BBG class of BzKs(uRJs) at 
z ~ 2{z ~ 1) and call these BBGs selected by including 
the AKARI NIR photometry as extended BBGs while the 
BBGs selected only with the ground-based photometry as 
classical BBGs. 



Fig. 50. The same as figure 48 for all the BBGs, includ- 
ing the extended BBGs. Top: For all the BzKs, includ- 
ing uzN3s and BzN3s. For BzKs-|-uzN3s-|-BzN3s also clas- 
sified as uRJs-|-uRN3s-|-BRN3s at 1.2 < z < 1.4, their sym- 
bols are changed to green triangles from red circles. For 
BzKs-|-uzN3s+BzN3s also classified as uVis at 0.4 < z < 0.8, 
their symbols are changed to blue stars from red circle. 
Middle: The same as on the left for all the uRJs, including 
uRN3s and BRN3s, as the Top. For uRJs-|-uRN3s-|-BRN3s 
also classified as BzKs-|-uzN3s-|-BzN3s at 1.4 < z < 2.5, their 
symbols are changed to red circles from green triangles. For 
uRJs-|-uRN3s-|-BRN3s also classified as uVis at 0.4 < z < 0.8, 
their symbols are changed to blue stars from green triangles. 
Bottom: The same as the Top for the uVis. For uVis also 
classified as uRJs-|-uRN3s-|-BRN3s at 1.2 < z < 1.8, their sym- 
bols arc changed to green triangles from blue stars. For uVis 
also classified as BzKs-|-uzN3s-f BzNSs at 1.8 < z < 2.5, their 
symbols arc changed to red circles from blue stars. 

Appendix 2. Summary of classifications 

Hereafter, we will summarize the photometric classifi- 
cation schemes introduced in the paper and show their 
outline as following tree diagrams, in which asterisk *, 
circle o, and star ★ at the front of a category represent 
a classification without using redshift information, with 
using redshift information, and with substituting photo- 
metric color conditions for redshift selections. 



32 



Hanami et al. 



A. 2.1. Optical classifications 

z'-detected galaxy 
o EBO 

o Rest-frame F-detected galaxy 

red sequence 

blue eloud 



BBG 



uVi 

*s-uVi 

*p-uVi 

*qp-uVi 

uRJ+(BRN3,uRN3) 

*s-uRJ 

*p-uRJ 

*qp-uRJ 



BzK+(BzN3,uzN3) 

*s-BzK 

*p-BzK 

>i=qp-BzK 



• z'-detected galaxies: As shown in subsection 3.1, 
we selected ~56000 of them, performed their pho- 
tometry as summarized in table 2, and estimated 
Zphot,M^,,AY, and SFRopt-corr as described in sec- 
tion 4 and appendix 3. 

• EBOs: The Extremely Blue Objects were selected 
with the criteria of equation (2), which are possi- 
bly AGNs as discussed in subsection 4.2 and ap- 
pendix 5.1. 

• Rest-frame ^-detected galaxy: We studied their 
stellar populations, which are subclassified as red se- 
quence and blue cloud in the CMD with criteria (5) 
as discussed in subsection 4.2. 

• BBGs: The Balmer Break Galaxies were introduced 
for redshift subclassifications of uVis, (e-)uRJs, and 
(e-)BzKs and distinction between star-forming (s- 
BBGs), passive (p-BBGs), and quasi-passive (qp- 
BBGs) populations with the criteria (A2), (A3), and 
(A4) in appendix 1. 

A. 2. 2. NIR classifications 



z'-detected galaxy 
I *N3R 



IRBG 



-*N23 bumper 
-*N34 bumper 
-*N4R 



N3N4R 

N3Rs: The N3 Red galaxies are selected from the z'- 
detected galaxies by the criterion (6) in subsection 5, 
which roughly correspond to populations at z > 0.4. 
IRBGs: The N3Rs were subclassified to the N23, 
N34 bumpers, and N4 Red galaxies (N4Rs) by the 
criterion (7) as InfraRed Bump Galaxies with 1.6 /J,m 
IR bump features from their stellar emissions. 



• N23, N34 bumpers, and N4Rs: They approximately 
correspond to the uVis. e-uRJs, and e-BzKs, respec- 
tively. 

• N3N4Rs: The N3N4 Red galaxies were classified by 
equation (12), which are possibly bright AGN can- 
didates as discussed in subsection 6.2. 

A. 2. 3. MIR classifications 

z'-detected galaxy 

*N3 Red galaxy(N3R) 
MbN3R 



7^s-MbN3R 

sb-MbN3R 



-s/a-MbN3R 



*agn-MbN3R 
MmN3R#:6.3 

*pP-MmN3R 

sb-pP-MmN3R 

non sb-pP-MmN3R 



★fP-MmN3R 

sb-fP-MmN3R 

non sb-fP-MmN3R 



7k-MmN4R 



★MfN3R 



MbNSRs: The MIR bright N3Rs were selected from 
the N3Rs with the detection in more than two MIR 
bands, as described in subsection 6.1. 
s- and agn-MbN3R: The MbN3Rs were subclassified 
to starbursts (s-MbN3Rs) and AGNs (agn-MbN3Rs) 
by the MIR color criteria (8), (9), and (10) at 2 = 
0.4 — 0.8,0.8 — 1.2, and > 1.2, respectively, as shown 
in figure 6. 

sb- and s/a-MbN3Rs: The s-MbN3Rs were sub- 
classified to StarBurst-dominant (sb-MbN3Rs) and 
Starburst/AGN mixturc(s/a-MbN3Rs) by their IR 
SED fittings. 

MmN3Rs: The MIR marginally-detected N3Rs were 
detected only in one or two MIR bands with SN3, 
which may be mimics of MbN3Rs with fainter MIR 
emission. They were subclassified the pP-, fP-, sb- 
and non sb-MmN3Rs as discussed in subsection 6.3. 
pP- and fP-MmN3Rs: The MmN3Rs were subclas- 
sified to the Possibly PAH emitting (pP-MmN3Rs) 
and the faint PAH (fP-MmN3Rs) populations as 
the PAH 7.7 um emission from the former is pos- 
sibly dominant in 511, L15, and L18 bands at z = 
0.4-0.65,0.65- 1.5, and > 1.5 while that from the 
latter is not. 

sb- and non sb-MmN3Rs: The StarBurst-dominant 
populations (sb-MmN3Rs) have MIR SEDs which 
are reproduced by the S&K model with a low AGN 
mixture rate < 50% while non sb-MmN3Rs do not. 
MfN3Rs: The MIR faint N3Rs were not detected 
in any MIR bands with SN3 as discussed in subsec- 
tion 6.4. 
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Dusty IRBGs/LIRGs at z = 0.4 - 2 
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A. 2. 4- Combining classifications 

We have sometimes introduced subgroups by combining 
these optical, NIR, and MIR classifications. For example, 
the N3Rs, classified by the criteria of the IR bump in 
the NIR photometry, can be also subclassificd into s-,p- 
, and qp-BBGs with the criteria of the Balmer break in 
the optical photometry, which were applied for the N3 
Red BBGs as the exclusively classified BBGs from the 
N3Rs as seen in subsection 5.2. Thus, the combination 
of classifications at different wavelengths can help us in 
quickly studying galaxies. 



Appendix 3. 
populations 



Photometric redshifts of various 



We obtained the spectra of 83 N3Rs, in which MIR 
bright sources are subclassificd into 35 sb-MbN3Rs, 9 
s/a-MbN3Rs, and 14 agn-MbN3Rs. For these spectro- 
scopic N3Rs and MbN3Rs samples, figures 51 and 52 
show the redshift discrepancy Az = Zphot — Zspec- The me- 
dian accuracies of {Az) / {1 + Zspec) are 0.040, 0.041, 0.039, 
0.042, 0.034, and 0.21 for all the N3Rs, relatively bright 
N3Rs (7V3 < 19.9), relatively faint N3Rs (7V3 > 19.9), sb- 
MbNSRs, s/a-MbN3Rs, and agn-MbN3Rs, respectively. 

We can see no essential difference in the median accu- 
racy between the relatively bright and faint N3Rs, which 
means that the brightness does not essentially affect the 
photometric redshift estimations. We can also see that 
sb-MbN3Rs, preselected for z > 0.4 with the IR bump 
criterion N2 — N3 > —0.3, show good agreement between 
Zgpec and Zphot , which means that the photometric redshift 
estimation is still robust even under heavy extinctions in 
these MbN3Rs as long as one excludes low redshift dusty 
starbursts with the IR bump detection. 

The agn-MbN3Rs show the largest discrepancy in all 
the subclasses in the N3Rs. It confirms that AGN is one 
major cause for inducing deviations in the photometric 
redshifts, which is reasonable since we used only stellar 
population synthesis modeled SEDs as BC03 without in- 
cluding AGN SEDs. Even though 14 spcctroscopically ob- 
served agn-MbN3Rs, as candidates harboring AGNs, in- 
clude 8 outliers, all the outliers are spcctroscopically iden- 
tified as BL AGNs and the remaining 6 non-outliers out 
of BL AGNs show that their median accuracy is compara- 
ble to that of the sb-MbN3Rs as (A2)/(l + z^pec) = 0.042. 
It suggests that their photometric redshifts can be ac- 
curately estimated from the SED fitting with their major 
stellar emission of the host galaxies even harboring AGNs. 

Thus, the photometric redshifts for N3Rs and MbN3Rs 
can be used to reconstruct their redshifts. 

Appendix 4. Redshift distributions of galELxies 

Figures 53, 54. and 55 show the photometric red- 
shift Zphot vs. rest-frame H absolute magnitudes and 
monochromatic luminosities vLjj{7.7 /im). Following the 
conversion from rest-frame H absolute magnitude to 
with equation (All) in appendix 7, we can see that galax- 




spec 




Fig. 51. The same as figure 1 for the N3 Red BBGs. The 
colors and symbols arc the same as figure 6. The magenta 
square means BL AGN classified with the spectroscopy. Top: 
For all the spectroscopic N3Rs. Bottom-Left: For reltaively 
bright N3Rs with A^S < 19.9. Bottom-Right: For relatively 
not-bright NSRs with 7V3 > 19.9. 

ies with > IO^^Mq are detected as the N3Rs at z > 1 
in figire 53. The standard nomenclature LIRG refers to 
a luminosity class of galaxies with L ~ lO^^^^^i© and 
is distinct from starbursts L ~ lO^^^^^L© and ULIRG 
> IQ^'^Lq. Following the conversion from vLuiJ .7 fiva) to 
TIR luminosity with equation (16) in subsection 7, we can 
see that most of the AKARI MIR-detected galaxies, classi- 
fied as the s- and s/a-MbN3Rs and sb-MmN3Rs, have typ- 
ically a TIR luminosity of IO^^^^^Lq, which corresponds 
to those of local LIRGs. 

We can see also that the trend of agn-MbN3Rs is dif- 
ferent from those of sb- and s/a-MbN3Rs. 
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Fig. 52. The same as figure 51 for the MbNSRs; the sb-, 
s/a-, and agn-MbN3Rs from left to right, which are subclas- 
sified on the basis of their MIR SED features as shown in 
subsection 6.1. 

Appendix 5. Models for SEDs and colors 

A. 5.1. Model tracks for extremely blue objects 

As remarked in subsection 4.1, all the spectroscopic 
BL AGNs in the Zphot outliers are also the EBOs with 
{U ~V) < ~0.25{Mv + 22.0) + 0.7, which suggests that 
most of the EBOs are possibly AGN candidates as long as 
the redshift estimation is unaffected even by SEDs with a 
mixture of stellar and AGN emissions. In order to confirm 
this assumption, we reconstruct the evolutionary tracks 
with the BC03 models on the CMD and the Magnitude 
vs. Stellar Mass diagram as shown in figure 56. Most of 
the observed star-forming galaxies suffer from extinctions 
Ay ~ 1 as discussed in subsection 7.3. Thus, most normal 
star- forming galaxies with Ay — 1, except young phase 
ones with age < 0.1 Gyr, cannot appear in the EBO re- 
gion of {U-V) < -0.25(A/y + 22.0) + 0.7 on the CMD. 
Even without any extinctions, only objects in an early star 
forming epoch < 1 Gyr can have the color of the EBOs. 
It is not realistic to ingnore extinctions. Indeed, we found 
evidence for existing dust even in the EBOs with MIR 
emission as reported in subsection 4.1. 

A. 5. 2. Model tracks for B aimer break galaxies 

In order to determine the criteria for selecting the 
Balmer break galaxies as shown in appendix 1.1, we stud- 
ied the model tracks on two-color diagrams of uRJ , uVi, 
BzN3, uzN3, BRN3, and uRN3 as shown in figures 57, 
58, 59, 60, 61, and 62, respectively. 

Top-left panels of all the figures represent the tracks 
of constant star-formation (CSF) for ages of 10^^, 
0.1,0.2,0.5,1, and 2 Gyr, various reddening with Calzetti's 
extinction law, and solar metalicity. We can see that 
galaxies in such a redshift range with ongoing star for- 
mation are indeed expected to lie in the left-upper side 
region B^iTZ > fB^iR of the two-color diagrams. The dura- 
tion of the star-formation (age) has little influence on the 
bluer color B — fi, while the redder color n — TZ increases 
with age due to the development of strong Balmer/4000 
A breaks falling in the /i-band around these selecting red- 
shift z ~ {z)b^iR.- Very young bursts with ages less than 
10 Myr without underlying older stellar populations would 
be located around the left-bottom region just below the 
line of BnTl= fBi_i,Tz- 

We will note that the reddening direction is approxi- 
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Fig. 53. Photometric redshift Zpjiot vs. rest-frame H abso- 
lute magnitudes for N3Rs. The colors and symbols are the 
same as figure 6. Upper-Loft: For all the N3Rs. Upper-Right: 
The same as the Upper-Left, for the IRBGs classified as s-B- 
BGs. Middle-Left: The same as the Upper-Left, for the 
IRBGs classified as p-BBGs. Middle-Right: The same as the 
Top-Left, for the AKARI MIR detected IRBGs classified as 
sb-MbN3Rs. Bottom-Loft: The same as the Top-Left, for the 
MIR detected IRBGs classified as s/a-MbN3Rs. Lower-Right: 
The same as the Upper-Left, for the MIR detected IRBGs 
classified as agn-MbN3Rs. 



mately parallel to the critical line of criteria B^TZ — fe^iTi 
since dust extinction is approximately represented as a 
power-law function of the wavelength and it is implying 
that the reddening is similar in the B — fi and fx — TZ colors 
for star-forming galaxies in the selected redshift range. It 
means that the generalized criteria BjiTZ > fs^iR. for se- 
lecting BBGs are also robust from their dust reddening. 

Top-right panels of all the figures represent the tracks 
of simple stellar population (SSP) models after an in- 
stantaneous burst, in the two-color diagrams, for ages of 
10~^,0.1,0.2,0.5, 1, and 2 Gyr, with no reddening, and so- 
lar metalicity. Their features are also similar to those in 
the BzKs as the tracks are similar to those of star-formings 
at young ages and begin to move to the region of p-BBGs 
at older ages of > IGyr. Even though some tracks for 
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Fig. 54. Photometric redstiift Zphot vs. rest-frame H abso- 
lute magnitudes for the MmNSRs. The colors and symbols 
are the same as figure 6. Left: For the sb-MmN3Rs. Right: 
For the non ab-MmN3Rs. 

intermediate ages of ^ 0.5 — 1 Gyr are outside of the s/p- 
BBG regions, the extinction tends to make them redder 
as those found in the s/p-BBG region. Even though we 
should be cautious when applying SSP models with no 
reddening to the analysis for real galaxies, the results sup- 
port that the generalized schemes for selecting s/p-BBGs 
might still be robust. 

The left-bottom panels of all the figures show that the 
Bj/lZ colors for galaxies with various formation redshifts 
are exponentially declining SFRs (r — 0.3 and 1 Gyr), with 
no reddening and solar metalicity. The color evolution of 
galaxies formed at high redshifts is such that most objects 
move directly from the star-forming galaxy region to the 
passive galaxy region without crossing the bluer regions 
populated by z < 1.4 objects. 

The right-bottom panel of all the figures shows that 
the colors of normal galaxies are derived from the CWW 
templates of E-Sbc-Scd-Irregular galaxies (Coleman et al. 
1980), which fall outside the color area of the BBGs for 
z > 0.8(z > 0.4), respectively. It means that the gener- 
alized two-color criteria can work to exclude the lower-z 
interlopers. 

A. 5. 3. Model Tracks for dusty starbursts and A GNs 

In figure 63, the thin blue curves, thick solid black 
curves, and black filled circles represent the redshifted 
tracks of dusty starburst SEDs of models by Siebenmorgen 
& Kriigel (2007), SEDs of Arp220, M82 and Mrk231, 
and a model torus SED used to fit the SED of a heav- 
ily obscured type 2 QSOs SWIRE J104409.95-F585224.8 
(Polletta et al. 2006), respectively. As shown in fig- 
ure 63, all the dusty starbursts of S&K models, Arp220, 
and M82 appear in the right area while the AGNs, includ- 
ing Mrk231 as an obscured AGN, appear in the left area. 
Thus, dusty starburst dominant LIRGs can be distin- 
guished from AGN dominant LIRGs on the MIR color di- 
agrams with a boundary represented as solid lines, which 
are introduced with equation (8) at z = 0.4 — 1.0, equa- 
tion (9) at z = 0.8 - 1.4, and equation (10) at 1.0 < z < 2.8 
in section 6. 
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Fig. 55. Photometric redshifts Zj^i^^t vs. monochromatic lu- 
minosities i/Li,{7.7 fim) at the rest frame 7.7 fim. The colors 
and symbols are the same as figure 6. 

A. 5.4. IR SED models 

With integrating a fitted model SED in rest-frame 8- 
1000 /.tm wavelength and correcting for the photometric 
redshift Zphot, we can obtain the total IR (TIR) luminosity 
Lju of the emission from the dusty star forming regions. 
There are various IR SED models. IR SED libraries can 
be prepared with three different kinds of schemes: 1) us- 
ing SED templates of known galaxies (Laurent et al. 2000; 
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Fig. 56. Left: The evolutionary tracks on the color-magni- 
tude diagram of My vs. {U — V) from the BC03 models. Red 
and black lines correspond to systems with M« = IO^^Mq 
and Mt = 10^^ Mq at 5 Gyr, repsectively, which are evolv- 
ing with exponential decaying SF with time scales of r =0.1, 
0.3, 1,2, 3,5, 15, 30 Gyr and Constant Star Formation (CSF) 
from the left to the right on the diagram. All of them are 
in the cases without extinctions. The open circles repre- 
sent ages of the systems of 0.1, 1, and 5 Gyr for all the 
models. An Arrow at the {Mv,U — V) = (—18,1) rep- 
resents color-magnitude correction with a mean extinction 
Ay ~ 1. Right: The same as on the left for the mass- 
magnitude diagram of My vs. M« from the BC03 models. 

Lutz et al. 2003; Spoon et al. 2004), 2) reproducing the 
SED with optical thin emission from the dusts assumed 
to be heated in a given radiation field without radiative 
transfer calculations (Chary & Elbaz 2001; Dale & Helou 
2002; Lagache et al. 2003), and 3) calculating radiative 
transfer to evaluate the emission reprocessed with dusty 
medium in a spheroidal shape (Silva et al. 1998; Efstathiou 
et al. 2000; Takagi et al. 2003; Dopita et al. 2006; Piovan 
et al. 2006b, a; Siebenmorgen & Kriigel 2007), or a general 
three dimensional shape (Kylafis & Bahcall 1987; Popescu 
et al. 2000), using ray tracing or Monte Carlo techniques 
(Bianchi et al. 2000). Even in scheme 3), the interac- 
tion between dust and radiation is treated consistently in 
Takagi et al. (2003) and Siebenmorgen & Kriigel (2007). 
Even though empirical conversions from a monochromatic 
IR luminosity have been sometimes used, their empirical 
laws have been frequently derived from the SED fitting 
(Bavouzet et al. 2008). When we choose to use them, 
we should take into account their conceptional differences 
and applicable situations along with their limitations. 

In order to extract information related to the physical 
properties in a dusty star-forming nucleus such as optical 
depth and dust density from the MIR SEDs with PAH 
emissions and Si absorption, in this paper, we have taken 
the S&K model. This can be used to prepare prepare 
MIR SEDs not only with various PAH emssion strengths 
but also Si absorption depth mainly parametrized with 
radiation field strength and optical depth, respectively. 
Individual SED features in the S&K model depend on 
the following parameters; total luminosity Lq, size 
of starbursting nucleus R kpc visual extinction Ay-sK, 
ratio of OB star luminosity with hot spots to the total 
luminosity Lqb/L^"*', and hydrogen number density in 
the hot spots n^^ . Even though the S&K model depends 
on the luminosity we have scaled luminosities and 
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Fig. 57. The evolutionary tracks in the uRJ diagram from 
theoretical models the same as the Fig. 8 for the BzK dia- 
gram in Daddi et al. (2004). The top-left panel shows con- 
tinuous star formation model tracks for ages from 1 Myr to 2 
Gyr and for E(B-V)=0, 0.3, 0.6. Top right panel has simple 
stellar population models for ages from 0.1 to 2 Gyr and no 
reddening. The bottom-left panel shows evolutionary mod- 
els with various formation redshifts and SFR timcscalcs and 
no reddening. At decreasing redshifts the tracks generally 
turn from bottom-left to top-right. The bottom-right panel 
show colors for the local templates of various galaxy types 
from Coleman et al. (1980). All models arc plotted for the 
range 0.8 < z < 1.8, except in the bottom-right panel where 
plots arc for < 2 < 0.8. The limits and color ranges of all 
four panels reproduce those of figure 45, for direct reference. 

used only their features in the MIR SED fittings. The 
parameters i?, and Ay-SK have some degeneracy for 
variations of the MIR SED features. As shown in the 
comparison between figures 63, 6, and 11, the mod- 
els with a hot spot parameter n'"' = lO^cm"-^ can cover 
typical MIR colors of MbN3Rs and MmN3Rs while the 
model with n''* = 10'*cm~'^ cannot. Thus, we have fixed 
the hot spot parameter with n'*'' = lO^cm"'^. And, we have 
used its 27 models fixed with R — 3kpc and parametrized 
with L*°* = 10i° i,10" \ and IO^^.Il^^ ^i/;Sif = 2.2,6.7, 
and 17.9, and Lqb/L*"* = 0.40,0.60, and 0.90, which are 
treated as only qualitative parameters to reproduce va- 
rieties in the MIR SEDs, not as quantitative parameters 
describing an exact physical state in the starburst nucleus. 

In fact, the conversion from vL^^ ^ to Ljr depends on 
the IR SED models used in the fitting. For an exam- 
ple, as shown in figure 64, for most of the sb-MbN3Rs 
and MmN3Rs, we can see that the estimated Ljfi-sK 
from the SED fitting with the S&K model is approxi- 
mately larger by a factor of three than those with the 
Chary & Elbaz model (hereafter C&E model). The IR 
emission at > 10 /im in most of the S&K model SEDs 
dominate more than those in the C&E model, which 
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Fig. 58. The uVi diagram the same as figure 57 except 
the redshift range. All models are plotted for the range 
0.4 < z < 1.2, except in the bottom-right panel where plots 
arc for < z < 0.5. The limits and color ranges of all 
four panels reproduce those of figure 46, for direct reference. 



Fig. 59. The BzN3 diagram the same as figure 57 except 
the redshift range. All models are plotted for the range 
1.4 < z < 2.5, except in the bottom-right panel where plots 
are for < z < 1.4. The limits and color ranges of all 
four panels reproduce those of figure 49, for direct reference. 



causes the TIR luminosity Lir-sk estimated from the 
S&K model to be larger than the Ljfi-cE from the C&E 
model. The S&K model tends to derive larger TIR lu- 
minosity ratio of Ljj^.sx/uL,^t y when SEDs are fitted 
well with larger Ay-sK even for samples with the same 
vL^T j. On the other hand, the C&E model derives al- 
most the same TIR luminosity ratio Liji-ce/vL^'j 'j for 
them. Thus, the difFercncc among the three model pa- 
rameters for Ay-sK appears as the split multi-sequences 
in figure 64. Figure 64 shows multi-sequences with differ- 
ent ratios Lir-sk / Lir;CE and Ay-^sK more clearly than 
figure 17. 

In order to reproduce the dependence of the optical 
depth to derive TIR luminosity in the paper, we have used 
the templates in the public library by Siebenmorgen & 
Kriigel (2007) while the analysis with Takagi et al. (2003), 
for MIR detected objects in the field, are also reported in 
other papers (Takagi et al. 2007, 2010). 

Appendix 6. Error estimations 

With photometric redshifts, we have derived physical 
values for the detected galaxies as absolute magnitudes M, 
monochromatic luminosities vL^, stellar mass M*, and 
SFR. The physical values include errors from not only 
their photometric fluctuation Am in the observed band 
but also uncertainty of their photometric redshift. In their 
error estimations, we have assumed that the former and 
the latter are independent of each other. 

Absolute magnitudes M is defined with luminosity L,j, 
observed magnitude rrii, and flux Fi, at the observed fre- 
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where is the luminosity distance. Thus, since 
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dL{zY /{I + z), we estimate the error of absolute magni- 
tudes M as 
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where we used d{\ogiQx) = (logg 10) ^dx/x = QAiMx/x 
and Am^ = 1.09AF^/i^^. 

For the monochromatic luminosities vL^ and the stellar 
mass M*, with the similar schemes, their errors can be 
derived as 
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where we used the error Arriif in K band for the stellar 
mass error estimation. 

In the error estimations for SFRjjv and SFRir, we 
have used errors Am^ in the u band and the IRC band 




Fig. 60. The uzN3 diagram the same as figure 59 except 
the redshift range. All models are plotted for the range 
1.4 < 2 < 2.5, except in the bottom-right panel where plots 
are for < 2 < 1.4. The limits and color ranges of all 
four panels reproduce those of figure 49, for direct reference. 



Fig. 61. The BRN3 diagram the same as figure 57 except 
the redshift range. All models are plotted for the range 
0.8 < 2 < 1.8, except in the bottom-right panel where plots 
are for < 2 < 0.8. The limits and color ranges of all 
four panels reproduce those of figure 49, for direct reference. 



related to the rest-frame 7.7 /im, respectively. 

Appendix 7. Stellar mass from Optical/NIR 
SEDs 

The stellar mass M» is estimated in a standard way with 
the optical- NIR SED fitting for all the z'-detected galaxies 
and an alternative way from the AKARI/NIR photometry 
for the IRBGs as a subclass in the z'-detected galaxies. 

The rest-frame H band wavelength corresponds to the 
rest-frame 1.6-yum IR bump where the emission is dom- 
inated by the stellar components in galaxies. It means 
that the absolute magnitude in rest-frame 77-band Mh 
represents the underlying stellar mass content M,,. This 
provided us with two ways to estimate the Mh with our 
dataset; 1) the direct estimation from the AKARI/IRC 
photometry in the rest-frame ff-band at 1.6 /xm, and 2) 
the SED fitting with the ground-based photometric data. 

First, the N2,N3, and NA photometry can directly de- 
tect the 1.6-/im bump from the AKARl/lRC detected 
IRBGs even at z > 0.4. Using the interpolation scheme 
with equation (4) for the N2,N3, and NA photometry 
with Aq = 1.6 /xm and z = Zphot, we can estimate the abso- 
lute magnitude MH-Ninizphot) corresponding to the rest- 
frame i7-band directly from the iV2,A3,andA4 magni- 
tudes of the AKARI/IRC detected IRBGs. 

Second, the model fitting of Bruzual & Chariot (2003), 
with the whole ground-based u*BVRi'z'JKs SED, de- 
rives the absolute magnitude Mn-^opt with our small mod- 
ified version of the hyperz as instantaneously deriving the 
photometric redshift Zphot in subsection 4.1. It can ap- 



ply for all the z'-detected galaxies. Thus, the derived 
MH-opt{zphot) from the optical-NIR SED fitting with the 
ground-based photometric data is an alternative to deter- 
mine the stellar mass M*. 

The absolute magnitude MH;opt{zphot) from the 
optical-NIR SED fitting is well correlated with that 
MH-NiR{zphot) from the N2,N3, and NA photometry as 
shown in figure 65. Note this rest-frame Mn^NiRizphot) 
from the AKARI/NIR data can be approximately inde- 
pendent to the MH-opt{zphof) from the ground-based op- 
tical/NIR data, except for the photometric redshift Zphot 
derived from the optical dataset. Figure 65 presents not 
only the consistency between the ground-based NIR and 
the AKARI/NIR observations, but also accuracy for the 
photometric redshift from the SED fitting only with the 
ground-based data. 

Weak age dependence between M^, and Mh can be seen 
in figure 66, in which the deep dark and faint symbols 
represent old and young populations in the red sequence 
and blue cloud), respectively, in the {U — V) — Auv as 
discussed in subsection 9.1. 

Even though the relation between an absolute magni- 
tude M in a band and the stellar mass depends on 
the Initial Mass Function (IMF) and the star formation 
histories, they can be approximately described as 

i°Kio^)--«-4(^'^-^^")' 

where M^^ is the rest-frame absolute magnitude cor- 
responding to the stellar mass of 10 ^^Mq. As shown 
in figure 66, for the Mh in the H band, we can take 
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Fig. 62. The uRN3 diagram the same as figure 59 except 
the redshift range. All models are plotted for the range 
0.8 < z < 1.8, except in the bottom-right panel where plots 
are for < ^ < 0.8. The limits and color ranges of all 
four panels reproduce those of figure 49, for direct reference. 



M^^ = M}} = -23.6. Even though we mainly used the stel- 
lar mass derived directly from the optical-NIR SED 
fittings in the following, the conversion from Mh to M* 
is useful as a quick estimation of stellar mass M*. 

We also checked the relation between absolute rest- 
frame V-band magnitude My and stellar mass M, [My- 
M» diagram) derived from the optical-NIR SED fitting as 
shown in figure 67. The EBOs and the BL AGNs appear 
in the right region in the My-M^ diagram, where normal 
star forming galaxies can appear only at the early phase 
< 1 Gyr. Their results in My-M* is consistent with their 
one on the CMD in subsection 4.2 as the EBOs are can- 
didates of AGNs (see also appendix 5). After excluding 
these EBOs, we obtained a relation between My and M*, 
approximated as equation (All) with taking My^ = —22, 
as shown in figure 67. This conversion from My to Af^ is 
useful in making mass-color diagrams in section 9. 
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Fig. 64. Comparison between TIR luminosities Ljij esti- 
mated from the C&E model and the S&K model for the 
sb-MbN3Rs and the sb-MmN3Rs. Symbols are the same as 
in figure 6 except their colors. Blue, green, and red symbols 
represent the redshifts of 0.4 < z < 0.8, 0.8 < z < 1.2, and 
z > 1.2, respectively. Large, medium, and small symbols in- 
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of ^V;Sif = 2.2, 6.7, and 17.9, respectively. Left: For the 
sb-MbN3Rs. Right: For the sb-MmN3Rs. 
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Fig. 65. Comparison of the rest-frame J?-band magnitudes 
for the N3 Red BBGs between M^f.jvjii and Mfj-^pt, which 
are derived from the AKARI NIR photometry and the 
ground-based optical-NIR SED fitting, respectively. Symbols 
are the same as in figure 6 except their colors. Blue, green, 
and red symbols represent the redshifts of 0.4 < z < 0.8, 
0.8 < 2: < 1.2, and z > 1.2, respectively. 
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Bottom-Right: For the agn-MbN3Rs. 
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